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ABSTRACT 


Boltzmann's  equation,  with  a  collision  term  differing  somewhat  from 
the  usual  and  taking  long-range  interactions  into  account,  is  shown  to  be  a 
correct  basis  for  an  analysis  of  physical  phenomena  in  (ideal)  ionised 
gases.  All  macroscopic  quantities,  including  the  electromagnetic  field, 
are  identified  with  statistical  averages  of  the  corresponding  microscopic 
quentities . 

We  then  consider  the  inherent  tendency  of  the  system  towards  sta¬ 
tistical  equilibrium,  i.e.  minimum  or  zero  entropy  production.  A  new 
proof  is  developed  for  the  H -theorem,  valid  when  the  intermolecular  forces 
are  mainly  of  the  long-range  type.  The  characteristic  times  for  the  diffe¬ 
rent  relaxation  processes  are  computed,  and  the  equilibrium  conditions  for 
a  rotating  gas  in  external  and  internal  force  fields  are  deduced  by  exact 
methods . 

The  general  theory  of  transport  phenomena  is  then  considered.  The 
transport  equations  are  derived,  and  it  is  shown  that  to  get  a  complete  set 
of  equations  it  is  necessary  to  postulate  "phenomenological  relations" 
between  "fluxes"  and  "forces",  these  being  defined  in  accordance  with  the 
conditions  of  equlibrium. 

Making  specific  assumptions  about  the  properties  of  symmetry  of  the 
gas  when  a  magnetic  field  is  present,  we  derive  the  proper  form  of  the  co¬ 
efficient  tensors  with  the  aid  of  the  theory  of  Cartesian  tensors .  The  usual 
assumption  that  the  tensors  are  isotropic  is  not  justified. 

Some  important  features  of  irreversible  processes  are  discussed.  We 
then  proceed  to  solve  Boltzmann's  equation  by  successive  approximations. 
The  result  thus  obtained  from  microscopic  theory  are  in  complete  accord 
with  the  results  obtained  earlier  by  macroscopic  reasoning. 
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INTRODUCTION 


The  purpose  of  the  present  paper  is  to  give  £  presentation  of  certain 
physical  problems  and  mathematical  methods  encountered  in  that  branch  of 
physics  which  deals  with  ionized  gases  in  electromagnetic  fields. 

Our  aim  is  the  description  of  non -equilibrium  states  of  the  system, 
that  is  phenomena  of  relaxation  and  transport.  Considering  the  macrosco¬ 
pic  behavior  of  the  system  to  be  identical  with  the  averaged  microscopic 
behavior,  we  shall  build  our  theory  on  the  most  general  microscopic  basis. 
This  is  not,  as  almost  invariably  presupposed,  the  Boltzmann  equation,  but 
the  more  fundamental  Liouville  equation.  From  the  latter  equation  one  may, 
however,  under  certain  conditions  deduce  the  Boltzmann  equation. 

The  microscopic  description  permits  us  to  draw  a  number  of  interesting 
conclusions  concerning  the  behavior  of  the  ionized  gas.  In  addition  it  per¬ 
mits  an  exact  deduction  of  the  macroscopic  transport  equations. 

Special  attention  must  be  given  to  the  inherent  tendency  of  the  system 
towards  statistical  equilibrium.  While  the  intermolecular  mechanism  of 
interaction  cannot  be  expected  to  have  any  influence  upon  the  form  of  the 
final  equilibrium  conditions,  it  is  essential  for  the  rate  of  approach  to  this 
state.  When  computing  the  relaxation  times  of  the  system,  full  account 
must  be  taken  of  the  long-range  intermolecular  forces.  Such  forces  are 
predominant  in  ionized  gases,  a  fact  that  has  not  always  been  given  the  de¬ 
served  attention. 

The  form  of  the  equilibrium  conditions  is  essential  to  the  definition  of 
the  "fluxes"  and  "forces"  that  are  linearly  connected  in  the  phenomeno¬ 
logical  relations.  As  we  shall  see,  the  existence  of  these  relations  is  ne¬ 
cessary  in  order  to  get  a  complete  set  of  equations.  .This  is,  consequently, 
in  general  only  possible  for  states  sufficiently  near  equilibrium.  The  form 
of  the  equilibrium  conditions  has  most  often  been  taken  over  from  macro¬ 
scopic,  thermodynamical  theory.  Instead  we  shall  deduce,  from  our 
microscopic  equations,  equilibrium  conditions  that  are  somewhat  more 
general  than  the  usual  ones. 

Many  authors,  when  dealing  with  ionized  gases,  simply  assume  the 
validity  of  the  usual  macroscopic  equations  of  classical  fluid  dynamics, 
and  forget  that  an  isotropic  medium  has  been  presupposed.  It  is,  however, 
an  important  fact  that  a  magnetic  field  will  introduce  an  anisotropy  in  & 
medium  composed  of  ionized  molecules,  giving  it  a  rotational  symmetry. 

There  are  thus  at  least  four  frequent  defects  in  earlier  presentations 
of  the  theory  of  ionized  gases,  namely: 

a)  The  founding  of  the  theories  on  Boltzmann's  equation,  which  should 
be  replaced  by  Liouville 's  equation. 

b)  The  failure  to  take  into  account  the  long-range  intermolecular 
forces . 

c)  The  too  special  equilibrium  conditions  assumed. 

d)  The  failure  to  take  into  account  the  anisotropy  of  the  medium 
introduced  by  a  magnetic  field. 

These  are  considered  to  constitute  a  sufficient  reason  to  give  this 
presentation  of  the  theory. 

Our  aim  will  then  be  to  remove  these  defects  and  build  up  a  consistent 
theory  from  a  very  general  microscopic  base. 


APPROXIMATIONS 


The  only  approximations  to  be  made  from  the  outset  are  the  following: 

a)  No  quantum -mechanical  effects  are  considered. 

b)  Relativistic  effects  are  neglected,  and  the  correlation  region  for  the 
molecules  are  assumed  to  be  sufficiently  small  so  that  the  inter¬ 
actions  inside  this  domain  may  be  considered  instantanous . 

Other  approximations  to  be  made  in  the  course  of  development  will  be 
explicitly  stated  at  the  appropriate  place  in  the  text.  They  are  in  main: 

c)  For  the  greater  part  the  treatment  will  be  restricted  to  "ideal  gases 
the  mean  kinetic  energy  being  much  greater  than  the  electrostatic 
energy  of  interaction. 

d)  The  state  of  the  system  will  usually  be  supposed  to  be  near  equili¬ 
brium.  That  is,  the  deviation  in  the  macroscopic  quantities  from 
the  equilibrium  conditions  of  Section  16  are  small.  This  validates 
a  first-order  theory. 

e)  The  effects  of  the  surroundings  of  the  system  will  usually  be  repre¬ 
sented  by  external  macroscopic  force  fields. 


NOTATION 


Our  notation  mainly  follows  that  of  Chapman  and  Cowling  [4]  .  The 
vector  notation  of  Gibbs  is  adopted.  Scalars,  vectors  and  tensors  of  the 
second  (and  occasionally  of  higher)  order  are  denoted  by  the  following  re¬ 
presentative  symbols: 

s  ,  v  ,  T  . 


▼ 

V  :  V 

the  del 


Simple  scalar  product  is  denoted  by  a  dot,  e.g.,  u*v  or 
V  .  Repeated  scalar  product  is  denoted  by  more  than  one  dot,  e.g. 


Instead  of  the  usual  symbol  7 
or  nab  la  operator,  where  r 


the  symbol 


d  r 


is  adopted  for 


is  the  position  vector. 


I  is  the  unit  (second  order)  tensor  with  components  4ik  (Kronecker 

delta).  The  trace  or  divergence  of  a  second  order  tensor  U  is  defi¬ 
ned  as 


=  U  :  It  (summation  over  i)  . 

o 

From  any  tensor  U  a  n  on  -  d  i  v  e  r  g  e  nt  tensor,  denoted  by  U  , 
may  be  derived: 


o 

U 


u 


I  . 


If  the  elements  of  a  tensor  U  are-  a  new  tensor  U  ,  the  con¬ 
jugate  to  U,  is  defined  having  the  elements  . 

Any  tensor  U  may  be  written  as  the  sum  of  a  symmetrical  and 
an  antisymmetrical  tensor: 


U  =  i(U  +  V)  +  j  (  V  -  V  )  . 


The  other  symbols  used  will  be  defined  in  the  course  of  development. 


PART  I  . 


GENERAL  THEORY 


1.  INTRODUCTORY  REMARKS 

When  we  undertake  any  given  physical  description,  the  ayatam  we  wish 
to  describe  must  be  set  off  from  the  surrounding  mass  of  all  physical  objects, 
i.e.  the  anviroiwant.  In  the  following,  our  system  is  a  certain  quantity  of 
gas ,  and  the  environment  may  consist  of  othet  gases ,  fluids  or  solid  forms  ' 
of  matter.  Fields  of  force  will  not  be  considered  to  be  independent  objects, 
but  rather  as  being  attributes  of  material  systems,  existing  because  of  the 
internal  structure  of  the  latter  or  their  motion. 

A  system  may  be  divided  into  an  arbitrary  number  of  component  systems, 
but  the  partition  cannot,  or  need  not,  be  indefinitely  continued.  We  eventually 
arrive  at  a  usually  very  large  number  of  alamantary  ayatama,  where  we  no 
longer  need  to  know  the  internal  structure  of  each  system.  We  are  only  inter- 
ested  in  the  externally  directed  effect  of  the  elementary  system  on  its  environ¬ 
ment.  From  this  point  of  view  the  interior  of  an  elementary  system  does  not 
belong  to  the  proper  system. 

We  thus  define  both  an  outer  and  an  inner  limit  to  our  system.  In  our 
case  the  elementary  systems  are  molecules ,  atoms  and  electrons ,  all  of 
them  possessing  a  certain  electric  charge  and  other  molecular  properties. 

The  atata  of  the  system  is  described  by  certain  coordinates  of  state, 
the  choice  of  which  depends  on  the  type  of  description  we  elect  for  our 
system.  We  are  often  content  with  a  macroscopic  description,  and  the  co¬ 
ordinates  then  are  quantities  that  refer  to  the  system  as  a  whole,  i.e.  pres¬ 
sure,  temperature,  volume,  etc. 

But  more  often  than  not  such  a  description  is  too  coarse,  and  we  must 
apply  a  microscopic  description.  Here  the  macroscopic  phenomena  consti- 
ture  nothing  but  the  sum  of  a  multitude  of  molecular  effects . 

The  most  detailed  description  possible  of  a  system  is  made  by  giving 
the  states  of  all  the  elementary  systems .  When  the  system  is  a  gas ,  this 
means  that  we  give  the  total  number  of  molecular  coordinates  of  state,  such 
as  position  and  velocity,  for  each  molecule  at  a  given  time.  Because  of  the 
large  number  of  elementary  systems  in  the  proper  system,  we  must  usually 
resort  to  a  statistical  description. 

If  we  are  to  describe  the  development  of  a  system,  we  must  know  how 
the  parts  of  the  system  interact  with  each  other  and  with  the  environment. 

We  must,  in  other  words,  know: 

a)  The  inwardly  directed  effects  of  the  environment  on  the  system 
(on  the  parts  of  the  system,  the  elementary  system). 

b)  The  outwardly  directed  effects  of  the  elementary  systems  on 
their  environments  (i.e.  the  other  elementary  systems). 

Concerning  the  latter  effects  we  must  postulate  that  the  internal  struc¬ 
ture  of  the  elementary  systems  is  unchangeable  and  independent  of  the  en¬ 
vironment.  If  this  is  not  the  case,  we  must  perform  a  further  subdivision 
into  new  elementary  systems.  The  externally  directed  effect  of  the  elemen¬ 
tary  systems  may  then  be  represented  by  constant  proper  fields  of  force  of 
different  types,  and  can  be  considered  as  a  molecular  property.  From  a 
macroscopic  point  of  view  this  will  take  the  form  of  different  types  of  fields 
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of  force  -  electrical,  gravitational,  etc.  Other  forces  are  of  the  magnetic 
type,  arising  as  a  result  of  the  electric  transport  of  charge  by  the  mole¬ 
cules  . 

We  may  expect  that  the  environment  can  be  divided  into  component  and 
elementary  fields,  in  the  same  way  as  for  the  proper  system.  From  this 
point  of  view  the  effects  a)  of  the  environment  on  the  system  are  the  sums 
of  all  the  separate  effects  of  the  elementary  systems  in  the  environment. 
These  effects  must  be  ascribed  to  the  proper  fields  of  the  elementary  systems 
and  the  forces  arising  from  their  motion. 

The  whole  problem  of  interactions  thus  reduces  to  that  of  mutual  effects 
between  elementary  systems  .  These  we  must  know  if  we  are  to  make  an 
exact  description  of  the  development  of  a  system.  This  really  means  that 
we  extend  the  system  so  as  to  include  part  of  the  environment.  The  more 
exact  we  wish  our  description  to  be,  the  more  we  must  extend  the  limits  of 
our  system,  inwards  and  outwards. 

By  making  certain  assumptions,  however,  we  can  often  describe  the 
physical  processes  of  a  system  without  detailed  knowledge  of  the  environ¬ 
ment.  We  may,  for  instance,  assume  that  the  system  is  sufficiently  far 
removed  from  all  environmental  influences,  so  that  interaction  may  be 
neglected.  If  this  is  not  the  case,  we  may  assume  that  the  effects  can  be 
represented  by  constant  or  time -dependent  macroscopic  fields  of  force,  by 
various  limiting  conditions,  and  so  on. 

We  now  proceed  to  the  treatment  of  gases,  particularly  ionized  gases. 


2.  MICROSCOPIC  AND  STATISTICAL  DESCRIPTION. 
BOLTZMANN'S  EQUATION 


Suppose  that  our  system  -  a  certain  quantity  of  gas  -  consists  of  several 
component  systems  or  gases,  which  we  denote  by  the  indices  a  =  1,2,  ... 
Each  component  system  is  in  turn  composed  of  a  large  number  of  identical 
elementary  systems  or  gas  molecules.  We  assume  that  the  internal  struc¬ 
ture  of  the  elementary  systems  is  constant,  and  that  each  molecule  is  sur¬ 
rounded  by  a  constant  proper  field  of  forces.  If  the  molecules  have  the 
electric  charge  ea  ,  the  potential  of  this  force  will  be  ea/r  . 

From  a  microscopic  point  of  view,  the  state  of  the  system  is  given  when 
the  position  r  k  and  the  velocity  c  ak  of  all  molecules  are  given  (k  = 

1,2,  ...  ,  v#  ,  where  va  is  the  total  number  of  molecules  of  the  type  a  ). 
Each  molecule  may  thus  be  represented  by  a  point  in  a  six -dimensional 
position -velocity  space,  the  (r.o)-space.  The  microscopic  density  of  these 
points  in  (r,  c)-space  is  given  by  [l]  , 

F#  (r,c,t)  =  Z  «(r-rak(t))  fi(e-eak(t)),  a  =  1,2 .  (1.1) 

where  6  is  Dirac's  delta -function. 

The  microscopic  density  in  position  space  is 

Na(r,t)  =  /Fa(r,c,t)do  =  Z  4(r-r  #k(t))  .  (1.2) 

Now  having  defined  the  microscopic  state  of  our  system,  we  are  inter¬ 
ested  in  finding  equations  that  describe  the  development  of  state  with  time. 
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The  functions  iJu  must  satisfy  the  equation  of  continuity  (Liouville's 
equation) , 


C 


(F 


= 

a  m 


(1.3) 


where  F-  (r.c.t)  is  the  total  microscopic  field  of  force.  This  is  com¬ 
posed  of  tne  external  field  T  _  ,  representing  the  effect  of  the  environ¬ 

ment  on  the  system,  and  another  force  term  which  represents  the  effects  of 
all  the  molecular  fields  of  force  on  a  particle  a  in  (r,o), 

P«m  =  F«o  +  *  *  <r'°'r(lk<,'>’  0  Pk<*'»  •  (1’4) 

Here  Fag  is  the  contribution  to  the  field  of  force  at  the  point  (r,c)  from  a 
particle  p  with  the  coordinates  (r.,  (t1),  e  ,  (t ')).  The  latter  term  can 
also  be  written  as  P  P 


r  j  /  f 
P  k 


ap 


(r.c.r'.c)  6(r'-r^k)5(c  -c'pk)dr'dc  = 


*  C-.c.r',e')Fp(r.,...t.)dr.dc.  .  (1.5) 

The  integration  extends  over  all  points  r'  with  the  exception  of  the  point 
r'  =  r.  which  must  be  excluded  as  we  cannot  include  the  contribution  from 
the  proper  field  of  the  particle  a  in  r  . 


From  this  microscopic  description  we  can  now  proceed  to  a  statistical 
description.  In  so  doing  we  regard  all  quantities  as  being  functions  of  an 
initial  configuration  of  particle  coordinates  (r#k(0) ,  «ak(0)  ,  k  =  1,2,  ... 
va  ,  a  =  1,2  ..  .)  ,  and  average  over  a  certain  set  of  such  configurations. 
We  denote  the  averaged  quantities  by  (  >  . 


# 


The  mean  quantity 
*a  (r.C.t) 


^Fa^  is  identical  to  the  familiar  distribution  function 
=  <F#>  =  (Fa(P,c,t)>  .  (1.6) 


We  may  interpret  fgdr  de  as  the  probable  number  of  particles  a  in  the 
element  of  volume  dr  dc  in  (r,  c)-space. 

In  a  similar  way  we  can  define  the  distribution  function  of  the  second 
order, 

fap(r,c,r\c',t,t')  =  <FaF^> 

faa(r,c,r\c',t,f)  =  <FaFa^>  -  6(r-r')6(o-C)  <F#). 

Here  dc  dr'dc'  is  proportional  to  the  probability  of  finding  a  molecule 

a  in  the  element  dr  dc  at  the  point  (r,c)  simultaneously  with  the  presence 
of  a  molecule  of  the  type  P  in  the  element  dr'dc'  at  the  point  (r',c')atthe 

retarded  time  t'  =  t  -  1  -  "  r-i  .  (We  assume  that  all  effects  of  force  propa- 

c 

gate  at  the  velocity  of  light  c.) 

The  definitions  of  the  functions  of  higher  order  are  closely  analogous  with 
expressions  (1 . 6)  and  (1.7). 


a  *  P 

(1.7) 

a  ■  0 
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In  the  lime  way  we  can  define  statistical  mean  values  of  all  quantities 
that  depend  on  the  particle  coordinates  T  .  and  C  .  .  If,  for  instance, 
we  average  the  equation  (1.3)  we  get 


df 

df 

df 

+  o  • 

dt 

TF+,«o 

a 

'  To 

+ 

do 

*  I*.? £=P  dr'<‘°'  *  0  ■  <>•«> 


In  Section  8  we  shall  show  that  the  last  term  of  this  equation  partly  gives 
origin  to  a  correction  to  the  external  macroscopic  field  of  force  F  ,  and 

d  f  ° 

partly  to  the  so-called  "collision  term"  — 1 -  °  Equation  (1.8)  may  there¬ 
fore  be  brought  onto  the  form 


c 


df 

a 

"3F 


+  F« 


a  =  1,2,  ...  (1.9) 


This  is  Boltanann'  a  aquation,  which  is  the  fundamental  point  of  departure  of 
almost  all  exact  treatments  of  the  physical  conditions  in  gases. 


3.  THE  ELECTROMAGNETIC  FIELD 


The  microscopic,  electric  and  magnetic  fields  are  composed  of  an  ex¬ 
ternal  and  an  internal  part, 


P  k 


P  pk’W  Pk', 


*  Bo(r>t)  +  H  V^'p^W' 


(1.10) 


Here  I  o  and  Bg  are  the  contributions  to  the  internal  fields  at  a  point  T 
from  a  molecule  of  the  type  P  with  position  r'pk  and  velocity  c'p  k  . 

They  may  be  derived  from  the  scalar  and  vector  potentials  4  a  ^  A 
defined  by  P  P 


e  a  efl 

♦  o(r.t)  =  - - 2 - S  -6 

*  Kr'Skl 


*8  *(r'-r'pk)6(C,-e'pk) 

- ^ - BiL  dr’  dc'  , 


r-r' 


(1.11) 


l  ea « Air  i  e.e'Mf'-i’flJMc'  ’C.J 
A.fr.t)  =  i  -£-£*-  =  i  /  -I - rl!s - dr' dc'  . 

*  I'-'jkl  I r  ”  r"l 

We  thus  have  the  electric  charge  and  current  densities  due  to  a  molecule  P : 
Qe  (r.t)  =  ep«(r-rpk)  =  /e  6(r-rpk)  4(e-c  pk)dc  . 

(1.12) 

Jep(».‘)  =  Vpk6<r-rpk>  S  /epc4(r-rpk)6<c-°pk>dc  • 

It  follows  from  the  formulas  (1.  11)  that  the  internal  fields  En  and  Ba 
defined  by  P  P 
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d 

Jr 


I  1  * 
c  dt  AP  ’ 


—  4>e=  ° 
or  p  c  dt  P 


Bp  =  Jr  x  Ap  ’ 


(1. 13) 


satisfy  a  set  of  equations  which  is  completely  equivalent  to  Maxwell's  equa¬ 
tions  for  macroscopic  fields  (Brittin  [2]  )  , 

7r  ’  *p  =  4*9ep  • 

31- 

as 

II 

O 

(1.14) 

Tr  x  Ep  =  mZ  dt"  BP  ’ 

—  v  X  =i-  —  K  +  —  i 
dr  *P  ”  c  dt  P  c  *  ep’ 

The  external  fields  XQ  and  B  ,  which  may  be  derived  from  the  exter¬ 
nal  potentials  A  and  A  ,  naturally  satisfy  the  corresponding  homogeneous 
equations.  It  follows  that  the  total  microscopic  fields  (1.  10)  satisfy  Max¬ 
well's  equations  with  the  respective  charge  and  current  densities 


'•““Si  "•»  =  t  J®p  Fp  (r',  ©',  t'Jdr'de' 


Jem  • 


(115) 


The  fields  X  and  S  may  then  be  derived  from  the  total  scalar  and 
vector  potentials  m 


♦m  *  +o  *  “  *  ♦»  +  Ff>  • 


es°' 

=  A  +  2  2  A«  s  A„  +2/ — E -  Ffl  (r*.  c',t')dr'dc'  . 

m  °  p  k  p  °  |r-r'  |  P 


(1.16) 


Instead  of  the  set  (1.  10)  we  can  also  write 


The  contribution  to  the  force  F  (1-4)  from  these  fields  is 

am'  ' 


(1.17) 
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e  „ 

P  =  - (1 

am  ma  ' 


m 


+  -  o  x  B  )  = 

c  rrr 


=  — [(*  +-cxB  )+  *  J(*A  +  I  e 

m_  L'oc  o'  fl*'P  c 


Bp)Fp  dr*  do1] 


(1.18) 


According  to  equations  (1.  17)  the  averaged  macroscopic,  electric  and  magne 
tic  fields  are 


dr'  do' 


(119) 


while  the  averaged  force  ^Fa  Fam\  which  is  included  in  equation  (1.8)  has 
the  form  ' 


<F« 


a  Fa  m) 


[£.(•, 


m 


+  i 

c 


c  x  B0)  +  2/(*p 


-  c  x  Bp  )fap  dr'  do'  ]  . 


(1.20) 


Hence  this  force  cannot  be  immediately  expressed  in  the  usual  way  by  the 
fields  X  and  B  . 

It  follows  from  the  above  equations  that  the  macroscopic  fields  E  and 
B  may  be  derived  from  the  potentials  <+m)  and  A  =  (A  \  ,  and 

satisfy  Maxwell's  equations  with  charge  density  <}  =  X  and^ 

J  »  (i  \  .  e  \  em/ 

’'e  \'em/ 


4.  THE  COLLISION  TERM 
6  i 

The  collision  term  — represents  the  effects  of  direct  interaction 

o  t 

between  molecules  and  their  proper  fields .  The  major  difficulties  in  solving 
Boltsmann's  equation  are  usually  connected  with  this  term. 

We  shall  now  discuss  the  conditions  for  a  gas  composed  of  two  kinds  of 
particles  with  positive  and  negative  electric  charge .  We  shall  see  in  Section 
4  that  the  correlation  length  for  interaction  between  identical  and  non- 
identical  particles  O  is  given  by 


when  e,  n  and  T  is  the  electric  charge,  density  and  temperature,  respec 
tively,  of  the  particles. 

Two  other  distances  of  interest  are 


Q 


o 


and  dsn 


(1.22) 
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Two  particles  that  collide  with  each  other  with  the  impact  parameter  Q 
will  on  the  average  be  deflected  by  an  angle  «/2,  and  d  is  the  mean 
distance  between  the  particles. 

If  we  now  assume  that  the  kinetic  energy  of  the  particles  is  much 
larger  than  the  mean  electrostatic  energy  of  interaction,  i.e. 


* 


2 

e 

$ 


Qq  «  d  , 


(1.23) 


it  follows  immediately  from  expression  (1.21)  that 

1  «  nD3  ,  d  «  D  .  (1-24) 

The  assumption  (1.23)  is  nearly  always  satisfied.  The  impact  parameter 
Q  for  deflection  at  a  right  angle  is  thus  much  smaller  than  the  mean  par¬ 
ticle  distance,  and  the  latter  is  much  smaller  than  the  correlation  distance 
D, 

qo  «  d  «  D  .  (1-25) 


Employing  the  usual  terminology  [3]  ,  we  let  the  term  "short-range 
collisions"  mean  interactions  between  particles  with  impact  parameter  in 
the  region  (0,  Qq).  As  Qq  «  d  ,  such  collisions  will  practically  always 

be  binary.  The  contribution  of  these  interactions  to  the  collision  term 

dc£« 

-  is  given  by  the  binary  collision  integral  [4]  , 

dt 

“-4r  •  ?(^yjf['«*'«V0V*^(c«VVKe,,d'>d*dv  <>-26> 

where  wfl ^  =  |(ca  -  Cp)  |  is  the  relative  velocity  between  two  particles 

a  and  0  before  collision.  The  quantities  ea  and  Op  are  the  velocities 

of  the  particles  after  collision,  and  are  functions  of  Oa,  Cg  and  b  .  The 

impact  parameter  b  -  the  distance  of  closest  approach  of  the  particles  i* 
the  absence  of  mutual  influences  -  and  the  polar  angle  c  characterise  the 
collision  geometrically .  The  region  of  integration  of  b  is  (0,  Qq  ) . 

When  the  impact  parameter  b  lies  in  the  region  ( Q0 ,  d)  we  no  longer 
may  assume  the  collisions  to  be  binary.  Encounters  with  more  than  two 
participating  particles  will  be  equally  probable. 

Interactions  with  impact  parameter  in  the  region  (d,  O)  may  no  longer 
be  characterized  as  collisions.  The  change  in  the  state  of  motion  of  a  par¬ 
ticle  arises  from  a  simultaneous  interaction  with  a  very  great  number  of 
particles .  In  view  of  the  large  number  of  particles ,  the  sum  of  the  many 
weak  effects  of  this  kind  most  often  becomes  much  larger  than  the  relative¬ 
ly  few,  but  strong  binary  collisions  (Section  13).  Such  multiple  interactions 
are  called  "long-range  interactions". 

We  shall  see  in  Section  10  that  such  multiple  interactions  can  be  repre¬ 
sented  by  a  "collision  integral"  of  the  Fokker -Planck  type,  first  derived  by 
Landau  [5]  , 
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d  f 
c  g 

d  t 


2*X 


£ 

P 


e«  *  0 


m, 


d 

7T 


/  <**«  fA  d£»  \ 

l\j.  JL  J. - —  -j— — )  do'.  (1.27) 

'  i  ffia  ic1  m#  do  /  '  ' 


Here  X  =  In  (d/ Q  q  )  is  a  cut-off  factor  which  arises  from  the  fact  that  we 
have  taken  into  consideration  interactions  with  impact  parameter  in  the 
region  (q„,  D).  We  have  thus  included  the  contribution  from  the  transition 
zone  (Q0,  d).  The  quantities  ea  and  e«  are  the  respective  charges  of 
the  particles  a  and  p,  while  ma  and  mp  are  their  corresponding  masses. 
Further,  fa  -  fa  (c)  and  fg  =  fp(c')  are  the  distribution  functions  for  the 
two  types  of  particles.  The  tensor  IT  is  defined  as 

U  =  -j  (u2  I  -  on)  ,  u  =  c  -  c 1  ,  u=|a|,  (1.28) 

u 

where  S  is  the  unit  tensor. 

Interactions  between  particles  with  impact  parameter  larger  than  the 
distance  D  may  be  characterized  as  plasma  oscillations  [6]  .  Their  effect 
must  be  included  in  the  macroscopic  force  term  in  Boltzmann's  equation. 

Before  turning  to  the  exact  derivation  of  Boltzmann's  equation  and  the 
determination  of  the  collision  term  (1.27),  we  shall  take  a  look  at  the  con¬ 
nection  between  the  macroscopic  and  the  (statistical)  microscopic  quanti¬ 
ties  of  a  gas . 


5.  MACROSCOPIC  AND  MICROSCOPIC  QUANTITIES 

The  distribution  function  fB  =  fa(r,e,t)  gives  the  density  of  an  aggre¬ 
gation  of  identical  molecules  a  at  the  point  (r,e)  in  phase  space  at  the  time 
t.  By  definition,  the  distribution  function  can  only  have  a  statistical  inter¬ 
pretation.  The  macroscopic,  measurable  quantities  in  which  we  are  inter¬ 
ested  may  be  identified  as  mean  values  over  velocity  space  of  microscopic 
quantities  with  the  weight  factor  fa  .  If,  therefore,  we  can  find  in  what 
way  fa  depends  on  r,  0  and  t  ,  we  have  the  exact  solution  to  the  rele¬ 
vant  physical  problem.  The  microscopic  quantities  mentioned  are  always 
various  kinds  of  molecular  properties,  such  as  number,  momentum ,  energy, 
mass  and  electric  charge.  We  let  ^  =  $(r,c,t)  denote  the  molecular  pro¬ 
perty  (generally  a  function  of  the  position  r  and  velocity  e  of  the  molecule, 
as  well  as  the  time  t),  and  define  the  mean  value  of  f  for  the  molecules  a 
in  the  following  manner, 

n«$«=  /4>afadc«  •  na  =  1  dca  •  (1-29> 

Here  na  is  the  density  of  the  molecules  a  in  coordinate  space.  In  general, 
both  na  and  $a  are  functions  of  r  and  t  ,  while  $  may  be  a  scalar,  a 
vector  or  a  tensor. 

If  the  gas  is  made  up  of  several  component  gases  (as  1,2,  ...),  we 
define  the  mean  value  $  =  $  (r,t)  for  the  whole  gas  by  the  expression 
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n  $  =  £  n«$a  =  J  J*afadoe  *  0-3°) 

where  n  =  Z  na  is  the  total  gas  density, 
a 

For  the  mean  velocity,  for  instance,  we  get 

no  =  I  naCa  =  I  /°afad°a  ,  (1.31) 

instead  of  which,  however,  general  usage  prefers  the  mats  Velocity  0o  > 
qco  =  ni5c  =  I  na5J7c^  =  J  <^ca  =  I  <^Cao  ,  (1.32) 

where  each  molecule  is  given  a  weight  which  is  proportional  to  the  mass  of 
the  molecule  ma  .  The  quantity  Qa  =  nama  is  the  partial  mass  density, 

and  Q  =  Z  Qa  is  the  total  mass  density.  The  momentum  of  the  gas  per 
ft 

unit  volume  is  evaluated  as  if  each  molecule  moved  with  the  velocity  cq  . 
The  mass  velocity  o  _  n  of  a  component  gas  coincides  with  the  mean  velo¬ 
city  . 

The  thermal  velocity  CB  of  a  molecule  in  a  composite  gas  is  general¬ 
ly  defined  in  relation  to  the  mass  velocity, 

Ca  =  0«  -  cQ  .  i.e.  I  qb  Cb  *  ®  .  (1-33) 

The  components  of  this  vector  are  denoted  by  (U,  V,  W). 

If  we  treat  each  component  gas  as  a  simple,  independent  gas,  it  is  more 
natural  to  define  the  thermal  velocity  as 

Ca  =  °a  -««o  '  i  e-  =  0  0-34) 

In  the  following  we  shall,  in  general,  use  the  first  definition,  but  we 
shall  also  have  occasion  to  make  use  of  the  latter  viewpoint. 

The  mean  thermal  energy  of  the  gas  (  $  =  JmC^)  is 

n(imC2)  =  Z  n  (imacj)  =  \  n  ,  *  =  kT  .  (1.35) 

This  is  to  be  considered  as  a  d*finttion  of  the  temperature  T  . 


6.  FLOW  OF  MOLECULAR  PROPERTIES 

We  shall  now  investigate  the  flow  of  molecules  through  a  small  "test 
surface"  with  area  dS  and  positive  normal  direction  a  [  4]  .  The  velo¬ 
city  of  the  surface  element  is  Of  .  The  relative  velocity  of  a  molecule 
with  respect  to  dS  is 

=  eH  -  e.  =  c  +  C  -  c, 

a  t  o  a  t  • 


(1-36) 
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We  consider  molecules  with  thermal  velocity  Ca  in  the  interval 
dCa  .  In  the  course  of  a  short  time  all  the  molecules  in  the  group  that 
were  originally  contained  in  the  volume 

dr  =  C r*'  •  dtdS  =  n  •  C r*1  dtdS 
a  a 

will  pass  through  dS.  The  number  of  molecules  thus  is 

fadCa  (a  •  C;*1  )  dtdS  .  (1.37) 

If  each  molecule  carries  along  the  molecular  property  $a(Ca)'  the 
total  amount  of  this  property  passing  through  the  surface  dS  during  the 
interval  of  time  dt  will  equal 

a  •  C"1  f  dC  dtdS  .  (1.38) 

a  ■  a  a  a 

The  total  amount  due  to  all  molecules  a  per  unit  time  and  unit  area  then 
becomes 


/  «*  •  C?1  <*afadCa  =  na($n  •  C?1  )  . 


(1.39) 


This  is  the  component  along  n  of  the  vector  (or  tensor)  na  (<J>a  Cat,!  )  , 
which  further  equals 


na(4>a  Co')  =  na^«Co)  +  na(co  ‘  cf}  *  O-40) 

The  last,  so-called  convective  part,  disappears  when  the  test  surface  dS 
moves  with  the  mean  velocity  of  the  molecules  c0.  In  this  case  we  obtain 
for  the  flux  of  the  property  , 

na  (4>«Ca  )  ,  (1.41) 

which  we  call  the  "flux  vector"  of  the  property  f  . 

The  flux  vector  for  the  number  of  particles  ($a  =  1)  thus  is 

n0C«  .  (1.42) 

and  disappears  for  a  simple  gas.  It  is  also  called  the  "diffusion  vector"  of 
the  molecules  a  . 

The  flux  vector  of  momentum  is  found  by  putting  $  =  mC  ,  which 
gives 

n  (macaca)  *  P0  (1.43) 

Another  term  for  this  tensor  is  the  "stress  tensor"  of  the  molecules  a  . 

In  the  same  way  the  flux  vector  of  energy  (  $  =  i  mC^)  is 


Ca)  =  la  •  (J-44) 

This  vector  is  usually  called  the  "heat  flow"  of  the  molecules  ci  . 
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The  expressions 

na0a  •  p  and  qa 

thus  give  the  flux  of  the  number  of  particles,  momentum  and  thermal  energy 
through  a  surface  which  follows  the  mean  motion  of  the  gas.  The  flux  of  a 
molecular  property  which  is  independent  of  the  velocity,  $a  =  Ka  ,  (e.g. 
mass  and  electric  charge)  is  evidently  given  by  the  vector 

n«KaCa  •  (1.45) 

The  total  contribution  for  all  component  gases  a  is  of  course  the  sum  of  the 
single  contributions, 

nKC  s  I  n^KgCg  , 


TP  =  nmC  C  =  I  nama®a®a  =  Z  jpft  ,  (1-46) 


q  =  n*mC2C  =  I  *  ma  c\  Cft  =  I  . 

According  to  equation  (1. 33)  the  mass  flux  vector  (Kg  =  ma  )  becomes 
zero,  while  the  electric  charge  flux  vector  (Ka  =  e#  )  becomes 

J*  =  neC  =  I  naeaCa  =  Z  naeaca  -  (  Z  na  ea  )c  .  (1-47) 

a  a  a  ° 

This,  therefore,  is  identical  to  the  conductive  current  density,  which  is  the 
difference  between  the  total  current  density  and  the  convective  flow. 

The  flux  of  momentum,  or  the  force  transported  in  the  gas  through  a 
surface  with  the  normal  direction  is 

Pi  =  *i-P  =  P’*i  =  nmCTc"  ,  (1.48) 

jp  being  symmetrical.  Part  of  this  force, 

2~ 

p.^  =  #.  •  p  -  «i  =  n(mC.  )  (no  summation)  ,  (1.49) 

works  perpendicularly  to  the  surface,  and  is  a  pure  pressure  force,  while 
the  remainder  works  parallel  to  the  surface,  and  is  a  pure  shear  force. 

For  instance, 

Pik  =  «i  •  P  *  «k  =  n(mCiC^)  ,  k  *  i  ,  (1-50) 

is  the  force  working  in  the  direction  on  a  surface  with  the  normal  . 

Writing  the  stress  tensor  on  matrix  form, 
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we  see  that  the  diagonal  terma  represent  the  pressure  forces,  while  the 
other  terms  represent  the  shear  forces . 

The  hydrostatic  pressure  is  defined  as  the  divergence  of  p  , 

P  =  y  (Pxx  +  Pyy  +  Pzz)  =  jPu  =  yP:I  = 

=  Inm  (U2  +  V2  +  W2)  =  jqC2  .  (1.52) 

In  the  case  of  equipartition  of  energy,  with  a  contribution  JkT  for  each 
degree  of  freedom,  we  get 

p  =  nkT  =  nl  .  (1.53) 

The  tensor  p  can  be  split  into  two  parts  (I  =  unit  tensor), 

°  1 

p  =  p  +  j  (p  :  I)  I  .  (1.54) 


The  first,  divergence -free  part  is  the  ahear^-atreaa  tensor  ,  and  the  other 
part,  the  pressure  tensor ,  represents  the  scalar  pressure. 


Analogously  with  equation  (1.51)  we  can  write  the  heat  flow  as  a  column 
matrix, 


Vu" 

r 

[*]  =  [nimC2cJ  = 

c2v 

= 

% 

C2W 

>q3. 

(1.55) 


where  the  terms  represent  the  energy  flux  in  the  three  coordinate  directions. 
We  may  use  a  similar  denotation  for  all  vectors,  e.g.  the  diffusion  vector. 

The  tensor  p  and  other  tensors  may  of  course  also  be  represented  by 
column  matrices,  where  the  nine  terms  follow  each  other  in  a  predetermined 
way.  We  shall  use  this  form  of  expression  on  a  later  occasion  (Section  21). 
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7.  THE  MACROSCOPIC  TRANSPORT  EQUATIONS 


If  we  perform  a  transformation  of  coordinates  from  the  set  (r,  e#,  t) 

then  Boltzmann's  equation  assu¬ 


to  the  set  (r, 
mes  the  form 

Ca  ,  t),  ca 

=  «a  -  «< 

df 

df 

do 

df 

-af  +Ca 

.  *L  +  (jp  . 
or  v  a 

—  )• 
dt  ’ 

a 

3T 

me 

dTT  Ca 


JL  e 

dr  o 


d  f 

c  a 

~3i~ 


(1.56) 


where 


d  .  d  .  d 
3T  "  dt  +  °o-ir  * 


(Two  points:  between  dyadics  means  that  scalar  multiplication  is  twice 
repeated. ) 


If  we  multiply  this  equation  by  $gdCa  where  $a(r,  C#  ,  t)  is  a  mole¬ 
cular  property,  and  next  integrate  over  velocity  space,  we  obtain  an  equation 
to  determine  the  course  of  development  of  the  corresponding  macroscopic 
quantity  n0$g, 
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(1.57) 


d$  .  d^n 
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The  right-hand  side  is  defined  by 


n 


« 


(1-58) 


and  gives  the  time  change  of  the  total  amount  of  properties  4>a  in  unit  volume 
due  to  the  collisions  experienced  by  the  molecular  carriers  a  of  the  property. 

The  other  effects  contributing  to  the  total  time  change  of  the  property 
in  unit  volume  are: 

a)  The  net  flux  of  molecules  out  of  or  into  the  volume. 

b)  The  way  in  which  $a  depends  on  position  and  velocity,  and  on 
time . 

If  we  substitute  C'a  for  Cc  from  equation  (1 . 34),  then  Boltzmann’s 
equation  (1.56)  and  the  macroscopic  equation  (1.57)  will  have  precisely  the 
same  form,  provided  that  we  replace  Cft  by  Cg  ,  and  «Q  by  C#0 
everywhere. 

We  can  now  insert  a  set  of  microscopic  quantities  in  equation  (1.57), 
obtaining  a  corresponding  set  of  macroscopic  transport  aquations.  The 
three  transport  equations  for  the  number  of  particles  ($  =  1),  momentum 
<♦  =  m  C)  and  energy  (  $  =  Jm  C2)  are  of  special  interest.  We  shall  re¬ 
turn  to  these  equations  later. 
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We  now  proceed  to  the  exact  derivation  of  Boltzmann's  equation  (1.9) 
and  the  determination  of  the  collision  term  (1.27). 


8.  DERIVATION  OF  BOLTZMANN'S  EQUATION 
The  averaged  Liouville's  equation  (1.8)  has  the  form 


if 


TT+°'  IF 


d£a  A 

♦r-  ■  -37 +  h  ■  dr'  *>'  -  0 


a  o 


(1-59) 


No  approximations  have  been  made  so  far.  In  the  case  of  only  electric  and 
magnetic  fields,  we  have,  according  to  (1.20) 


a  o 


=  i  c  x  BJ, 


ma  '  o  c 


*  sr  <*e  +  z  •  «  B» )  • 


As  a  result  of  interactions,  the  particles  will  not  move  independent  of 
each  other,  certain  correlations  will  exist  between  the  motions  of  the  in¬ 
dividual  particles.  This  physical  condition  may  be  mathematically  formu¬ 
lated  by  letting  fap  have  the  form  , 

fap(r,e,r',c',t,t')  =  4  (r.o.t)fp  (r'.c’.f)  +  ?ap(r,c,r'.c',t,t')  .  (1.60) 


The  correlation  factor  ?ap  is  only  important  in  a  region,  the  dimensions 
of  which  is  the  correlation  length  ( |r  -  I*'  |  5  D). 


The  distribution  functions  of  higher  orders  have  corresponding  forms. 
We  can  for  instance  write  f  apy  '• 

fapr  (r.e,r\c',r".c",t.t*,t")  =  fa  (r,c,t)fp(r',c'.f)fy(r",c",t'')  + 

+  f#(r,e.t)  ?pr  (r'.c'.T*'  ,c",t',t")  +  fp(r',c',t')  *ya  (r",o",r,c.t",t)  + 

+  fy(r",c".t")fap(p,«,r',c',t,t’)  +  fapy  (r,c,r',c',r",c",t,t',t")  .  (1.61) 
When  equation  (1.60)  is  substituted,  equation  (1.59)  assumes  the  form 

df-  df-  A 

75  +  h  ■  *  / t«p  *>'  1  o  •  (>•“) 


a  ,  o._ 

sr  +  °  •  tf  +  F« 


The  total  macroscopic  field  of  force  F  is  expressed  in  the  usual  manner 
by  the  averaged  electric  and  magnetic  fields  ■  =  ^  K  \  and  B  =  \ 

defined  by  equation  (1.  19)  (plus  forces  of  different  origin).  ^ 

We  permit  ourselves  to  assume  that  the  correlation  region  where  ^af) 
differs  from  zero  is  so  small  that  no  noticable  error  is  committed  by 
neglecting  the  time  retardation  and  other  relativistic  effects  within  this 
area.  We  may  then  neglect  the  magnetic  interaction  between  the  particles, 
and  only  take  account  of  the  electrostatic  interaction.  Equation  (1.62)  then 
reduces  to 


df« 
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To  to  p  ma  * 
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fapd»dc 


0  , 

(1.63) 
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where 


=  (C  +  7-  e  x  B) 


m 


Equation  (1.63)  is  the  starting  point  of  at  least  two  important  recent 
investigations,  the  first  by  Kadomtsev  [7]  and  a  later  one  by  Green  [8]  . 
Green  shows  how  one  can  deduce  macroscopic  equations  directly  from 
(1.63)  (by  taking  moments  in  velocity  space).  Kadomtsev’s  procedure  is 
not  quite  as  general  as  Green's,  but  it  gives  a  better  insight  into  the  micro¬ 
scopic  behavior  of  the  system.  In  Sections  8-10  we  shall  therefore  make  a 
further  study  of  the  implications  of  equation  (1.63)  along  the  general  lines 
of  presentation  of  Kadomtsev,  showing  how  the  usual  Boltzmann  equation 
can  be  deduced  from  equation  (1. 63)  (with  a  collision  term  differing  some¬ 
what  from  the  usual  ones).  Unlike  Kadomtsev,  however,  we  shall  allow 
magnetic  fields  that  vary  in  time,  our  equilibrium  conditions  (1.57)  are 
more  general  than  his  ,  and  our  results  are  in  part  valid  for  a  less  restric¬ 
ted  system. 

Later  on,  when  we  deduce  the  macroscopic  equations  from  Boltzmann's 
equation  we  shall  occasionally  have  the  opportunity  to  make  comparisons 
with  the  results  of  Green. 

Our  system  is  now  a  binary  mixture  of  electrons  (index  1)  and  positive 
ions  (index  2),  and  we  suppose  that  the  assumption  (1.23),  stating  that  the 
electrostatic  interaction  energy  of  the  particles  is  much  smaller  than  their 
kinetic  energy,  is  correct, 

e2/d  «  *  or  1  «  nD3  .  (1.64) 


In  this  case  all  correlation  factors  ?  are  small  quantities,  while  the  triple 
correlation  faPy  is  smaller  than  the  binary  factor  fap  by  one  order  of 
magnitude,  and  may  be  neglected. 


The  problem  is  now  to  find  fap  .  We  can  derive  the  equation  to  deter- 
e  fap  in  the  following  manner: 

We  write  down  Liouville's  equation  (1.3)  for  Fa  (r.e.t)  and  for 
Fa  (r'.St't').  We  multiply  the  first  equation  by  F'a  and  the  second  by  Fa 

.  5  J _ ] _  T _  A.X  •  _ _ _  -Li.  _  •  _ _ _ ./lai _ / _ r  _ _ 
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dt 


In  this  manner  we  obtain  an  equation  for  which  con- 
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dr"  do"  =  0 


(1.65) 


C  j 

Fflp  =  j^-(Ep  +  £  0  x  Bp)  is,  as  before,  the  force  per  unit  mass  on  a 
particle  a  at  the  "point"  (r,e)  due  to  a  particle  p  at  the  "point"  (r',0'). 
In  this  equation  we  substitute  the  expressions  (1.60)  and  ( 1 . 6 1 )  for 


f«p 


fQpy  and,  using  equation  (1.62),  obtain: 
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(1.66) 
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dir  y  P7  ®y 

Equation  (1.66)  is  very  complicated  and  must  be  solved  by  methods  of 
approximation.  We  first  require  that  the  equation  be  satisfied  in  the  station 

ary  case  (-^-  =  0)  by  substitution  of  the  quantities 
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I  • 


(1.67) 


We  next  require  that  in  this  case  the  equilibrium  conditions  (2. 91),  to  be 

deduced  later,  are  valid.  These  equilibrium  conditions  differ  from  those 

adopted  by  Kadomtsev  by  allowing  spatial  gradients  and  a  mean  velocity 

field  o  ♦  0  . 
o 

Our  next  step  is  a  transformation  of  variables  in  equation  (1.66)  from 
the  set  (r,c,  r' ,c' , t)  to  the  set  (r, C ,  r',  C,t),  where  C  =  , 

C'  =  c'  -  c'Q  .  On  inserting  the  expressions  (1.67)  into  equation  (1.66)  we 
obtain,  after  some  simplifications,  an  expression  which  is  a  polynomial  of 
the  second  degree  in  C  and  C.'  The  constant  term  and  the  coefficients  of 
the  second-order  terms  disappear  when  the  equilibrium  conditions  (2.91) 
are  introduced.  This  is  also  the  case  for  a  part  of  the  first-order  coeffici¬ 
ents.  The  remaining  part  must  also  disappear,  and  we  are  left  with  the 
equation 
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(1.68) 


On  putting  -e^  =  e£  =  e  »  n2  =  n2  =  n  *  Kadomtsev's  equation  results.  This 

equation  is  therefore  valid  for  more  general  equilibrium  conditions  than 
those  adopted  by  Kadomtsev. 
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Putting  in  equation  ( 1 . 68)  r  -  r1 
takes  the  form 
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(1.69) 


For  distances  Q  greater  than  the  distance  qq  -  e  V A  corresponding  to  a 

mean  deflection  */2  of  the  particles  when  they  collide,  we  can  neglect 
the  first  term  on  the  right-hand  side  of  equation  (1.69).  The  resulting  equa¬ 
tion  is  easily  integrated,  giving 
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(1.70) 


We  have  shown  that  the  arguments  of  Kadomtsev  may  so  far  be  extended 
to  a  more  general  situation.  As  the  further  computations  are  lacking,  we 
shall  nevertheless  be  satisfied  only  to  give  the  results  of  Kadomtsev.  We 
must  expect,  however,  to  get  similar  results  even  in  the  more  general  case. 

Assuming  with  Kadomtsev  that  - e ^  =  e 7  =  e ,  nl  =  n2  =  n’  (elec¬ 
trical  neutrality)  the  solution  of  equation  (  1.70)  is  givenby 


X«P«>  = 


ege(3 
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I  e“P/D 

Q 


where  D  is  the  Debye  shielding  distance, 


D  = 


(1.71) 


(1.72) 


and  at  distances  greater  than  D  the  correlation  between  the  particles  will 
rapidly  approach  zero.  The  condition  ( 1 . 64),  qD  =  e/i  «  d  ,  means 

-  2 

that  the  region  where  the  solution  (1.71)  is  not  valid,  q  <  ,  constitutes 

a  very  small  part  of  the  total  correlation  region.  We  shall  therefore  not 
worry  about  the  form  of  the  solution  within  this  region,  but  simply  neglect 
the  effect  of  collisions  with  impact  parameter  less  than  qQ  .  This  is  equi¬ 
valent  to  neglecting  the  short-range  collisions  (Section  4),  and  entails  that 
the  first  two  terms  on  the  right-hand  side  of  equation  (1.60)  may  be  dropped. 

We  shall  further  assume  that  the  deviations  in  the  state  from  the  equili¬ 
brium  conditions  (1.67)  are  slight.  We  may  then  insert  the  expression  fag 
(1.71)  into  the  right-hand  side  of  equation  (1.66)  which  is  small  in  terms  of 
l/nD*,  and  also  into  the  terms  on  the  left-hand  side  involving  the  macro¬ 
scopic  forces.  We  write  these  terms  as  the  sum  of  a  part  symmetrical  and 
a  part  antisymmetric al  in  q  =  r-r',  and  correspondingly  split  fga  into 
two  parts 

+  •  (1-73) 

are  then  defined  by 
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where  U(q)  =  -  e^ / Q  e  , 


(1.75) 


Q  =  *  -  r' 


What  is  now  left  is  to  solve  the  equations  (1.74)  and  (1 .75)  and  insert 
the  resulting  expression  (1.73)  into  equation  (1.63).  The  form  of  this  equa 
tion  is  then  determined. 

The  contribution  to  equation  (1.63)  from  the  symmetrical  part  f 
may  be  considered  as  a  correction  to  the  macroscopic  force  field.  This 
correction  has  the  form 
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and  is  in  general  a  function  of  r,  O  and  t  . 
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In  a  similar  way  the  antisymmetric  part  f '  '  will  give  a  contribution 
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to  the  equation.  This  contribution  we  define  as  th  collision  term. 

It  is  worth  noting  that  both  the  correction  AF  and  the  collision  term 
d  f  ° 

— —  will  vanish  when  the  correlation  between  the  particles,  represented  by 

<Pa  p1  vani*h«»- 


9.  THE  CORRECTION  TO  THE  MACROSCOPIC  FIELD 

We  shall  not  carry  out  the  detailed  solution  of  equation  (1.74)  as  has 
been  done  by  Kadomtsev.  His  results  may  be  summarized  in  the  following 
remark. 

The  existence  of  the  correction  term  AFa  may  be  explained  physically 
by  the  formation  around  each  gas  particle  of  a  correlated  "cloud”  of  particles 
of  opposite  charge  and  dimensions  ~  D3.  In  the  case  I  =  B  =  0  the  par¬ 
ticle  cloud  will  be  spherically  symmetrical,  but  when  X  ♦  0  the  "polarisa¬ 
tion"  AFa  will  make  itself  noticed.  When  the  magnetic  field  is  weak,  the 
polarization  will  be  along  -I,  while  in  the  opposite  case  the  displacement 
of  the  correlated  cloud  and  the  polarization  will  be  along  the  strong  magnetic 
field.  The  gas  thus  behaves  like  a  dielectric  medium. 

The  "polarization"  AFa  is,  however,  proportional  to  l/nD3  and  maybe 
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neglected  when  terms  of  this  order  of  magnitude  are  unimportant.  For  real 
gases,  on  the  other  hand,  where  the  interaction  energy  is  not  negligible,  the 
"virial  term"  plays  an  important  role. 


10.  THE  COLLISION  TERM 

{  \  d  f( 

We  now  try  to  solve  equation  (1.75)  for  9^(3  >  *n  order  to  find  —  y 

by  substitution  in  equation  (1.77). 

We  assume  that  the  distribution  functions  f<x  =  fa(l*.  o,t)  do  not  vary 
appreciably  in  the  correlation  region  (~D^),  and  moreover  that  they  are 
approximatively  constant  during  the  time  used  by  a  particle  to  cross  this 
region  (on  the  average),  we  may  neglect  the  dependence  in  time  and  position 
for  fa  ,  i.e.  fa  =  f a (c) •  We  further  assume  the  following  form  for 

J-J 
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Equation  (1.75)  then  takes  the  form 
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The  next  step  is  to  introduce  cylindrical  coordinates,  so  that  the  coordi¬ 
nates  for  q  are  (r,  f  ,  z).  The  z-axis  is  chosen  along  the  vector  o  -o'. 
We  then  get 
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which  can  be  immediately  integrated,  giving 
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We  now  substitute 


We  have  then  assumed  that  f  ^  ~  0  for  z 
9^  in  equation  (1.77),  getting  the  final  expression  for  the  collision  term, 
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This  is  the  collision  term  on  Landau's  form,  given  in  equation  (1.27).  This 
term  has  been  the  subject  of  some  controversies  in  scientific  literature. 
Enoch  [9]  has  recently  proved  the  equivalence  of  the  Landau  collision  term 
to  the  Fokker -Planck  term  on  the  assumption  that  binary  correlations  domi¬ 
nate. 
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With  this  collision  term  equation  (1.63)  takes  the  form 
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In  the  case  1  «  nD^  we  neglect  A  Pff  ,  and  equation  (1.83)  gives  us 
Boltzmann's  kinetic  equation  (1.9). 
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PART  II.  RELAXATION  PHENOMENA 
AND  EQUILIBRIUM  STATE 

1 1 .  THE  H -THEOREM  AND  THE  TENDENCY  TOWARDS 
STATISTICAL  EQUILIBRIUM 


When  we  coneider  the  statistical  state  of  a  system,  the  molecular 
interactions  are  of  decisive  importance.  Their  significance  has  been 
clearly  stated  by  Fowler  [10]  : 

"Any  possible  mechanism  (of  interaction)  left  to  act  by  itself  must  set 
up  and  preserve  the  laws  of  statistical  equilibrium." 

The  so-called  H -theorem  for  gases  was  first  proposed  by  Boltzmann 
in  1872  [11]  .  It  stated  that,  when  left  to  itself  without  interaction  with 
the  environment,  a  gas  would  approach  statistical  equilibrium  represented 
by  a  Maxwellian  distribution  function,  as  a  result  of  binary  collisions  among 
the  molecules  (E  is  tne  sum  of  kinetic  and  potential  energy,  E  =  ^  mc^  +  mf), 

f  =  n(£A)#/>  e_PE  .  P  =  1/kT  .  (2.1) 

This  distribution  was  earlier  found  by  Maxwell  in  I860,  but  at  the  time  his 
proof  was  not  considered  satisfactory. 

The  point  of  departure  for  Boltzmann's  proof  of  the  H -theorem  and  the 
subsequent  improvement  on  the  theorem  was  Boltzmann's  equation. 

In  the  Maxwell -Boltzmann  statistics  the  thermodynamic  "probability" 
for  a  system  is  given  by 
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(2.2) 


where  N.  is  the  population  number  and  the  statistical  "weight"  of  the 
i-th  phase  cell  in  the  six-dimensional  phase  space,  while  Ns  2  N.  is 

i  1 

the  total  number  of  molecules.  We  assume  W  to  have  an  increasing  ten¬ 
dency,  and  W  =  Waax  or  6W  =  0  characterize  statistical  equilibrium. 
The  latter  quantities  are  distinguished  by  Boltzmann's  law  of  distribution. 


P  =  1/kT  , 


(2.3) 


where  E^  is  the  energy  level  of  the  i-th  cell,  and  the  "sum  of  state" 

-PEt 

Z  =  I  G.  e 
i  1 

In  this  state  the  thermodynamic  expression  for  the  entropy  is  given  by 
S  =  In  W  (k  is  Boltzmann's  constant). 
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scopic  quantities.  What  interests  us  is  to  find  the  mtoroaoopio  juattfloatlon 
of  those  assumptions .  We  have  so  far  merely  assumed  that  the  thermodyna¬ 
mic  probability  of  a  system  has  a  growing  tendency,  with  the  result  that  our 
system  approaches  the  state  of  equilibrium.  With  the  aid  of  Boltzmann's 
kinetic  equation  we  can  now  directly  show  that  because  of  the  interaction 
between  the  molecular  systems  the  function  H  will  decrease  towards  a 
lower,  negative  limit  which  corresponds  to  thermodynamic  and  statistical 
equilibrium.  (The  existence  of  this  lower  limit  has  been  proved  by  Chapman 
and  Cowling  [4]  .)  The  state  is  then  stationary  with  a  Maxwellian  distribu¬ 
tion  function,  and  the  entropy  Sux  =  -kHain  corresponds  to  the  classic 

thermodynamic  entropy.  This  is  the  substance  of  Boltzmann's  H -theorem. 


12.  PROOF  OF  THE  H-THEOREM 


According  to  the  definition  (2.9)  the  time  variation  of  Ha  is  given  by 
the  expression 


dH« 

<Jt 


/  (1  +  lnfa) 


“f« 

IT 


de 


(2.10) 


df 

in  which  we  substitute  for  — from  equation  (1.83)  and  rearrange,  getting 


OH  d  H  „ 

tt  = ^0f«  lnfad0  +  /^  p;  )f«lnf«dc  *  (211) 


or 


*Ha  _  *cH«  d 

dt  "  dt  dr 


-  —  .  (n.c  In  f#)  +  n  (—  .  )  In  f 

t  dr  '  a  a  • 


de 


(2.  12) 


dcHa 

Here  is  the  change  in  Ha  because  of  the  molecular  interactions 

d  f 

c  « 


( J  dc  -  0  ,  see  Appendix), 


d_f 


d  In  f 

c  a  c  s 


=  /  (1  +  In  f#  )  -jp-  dc  =  na 


dt 


(2-13) 


If  we  write 
dH 


d  t  +  ’JJT  '  (n«  0  ln  )  =  Tt“  +  n«  ^  •  *s  )ln  fa  * 


(2.14) 


this  may  be  interpreted  as  an  equation  of  continuity  for  the  "entropy  density" 
Ha  .  This  equation  also  follows  from  equation  (1.57)  by  setting  fa  =  ln  fa  . 
The  vector  na  c  ln  fc  =  /  cf.lnf.de  is  the  "entropy  flow  density",  and 
the  right-hand  side  of  equation  (2.  14)  is  the  "entropy  source  density". 

As  shown  elsewhere  (Section  8)  the  effective  macroscopic  field  F£^ 
which  includes  AFa  may  be  a  function  of  c  .  We  can  therefore  not  imme - 
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There  ie  no  thermodynamic  definition  of  entropy  for  non-uniform  and 
non-stationary  states,  but  it  is  convenient  to  ohoose  the  equation 


S  =  k  In  W 


(2.4) 


as  an  equation  of  definition  even  in  this  case,  W  being  defined  for  all 
states.  This  procedure  seems  natural  as  both  W  and  S  are  assumed  to 
have  increasing  tendency,  and  because  the  entropy  for  thermodynamic 
systems  is  additive,  while  the  thermodynamic  "probability"  is  multiplica¬ 
tive.  (This  argument  has  been  criticized  by  Fowler  [10]  .) 

If  we  use  the  expression  (2.4)  for  W  and  Stirling's  formula 
In  N  !  sNlnN-N, 
we  can  easily  deduce  that 


/Ni  Ni\ 

In  W  *  -  Z  y  In  ^  j  Gj  +  N  In  N  . 


(2.5) 


The  statistical  "weight"  G.  is  by  definition  proportional  to  the  volume  of  ■ 
the  i-th  phase  cell  A^j  .  If  we  choose  the  proportionality  factor  equal  to 
unity  and  let  -*  h3  (h  is  Planck's  constant),  the  average  density  of 

particles  in  the  phase  cell. 


AN 


i 


Ni 

approaches  a  limit  f  ,  and  the  sum  (2.5)  can  be  replaced  by  an  integral. 

In  W  *  -  J  f  In  fd  +  +  Const. ,  d  +  =  dr  do  .  (2.6) 

Our  expression  (2.4)  for  the  generalized  entropy  thus  becomes 

S  =  S(t)  =  -kj  f  In  fd  ♦  + Const.  .  (2.7) 

If  we  have  a  composite  gas,  we  define  the  entropy  as  the  sum  of  the 
entropies  of  the  different  components  of  the  gas  (*  =  1,2,  .  .  .  ), 


S  =  2  S. 


=  -k  Jf#  In  f(  d 


(2.8) 


This  corresponds  to  the  total  thermodynamic  probability  being  equal  to  the 
product  of  all  component  probabilities . 

Instead  of  S, 
defined  by 


.a  it  is  customary  to  work  with  the  functions  H  a  Q  of  H#  , 


Hao(t)  *  /H«dr' 


H,(r.t)  =  /  f,l»fadc  , 


(2.9) 


i.e.  S. 


•kH 


ao 


Thus  H  is  proportional  to  the  negative  entropy. 

The  distribution  (2.  1)  now  correspond*  to  tae  equilibrium  distribution 
(2.3).  This  result,  however,  rests  on  several  assumptions  about  macro- 
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diately  put  —  •  Ftt  -  0  ,  even  though  the  contribution  from  the  Lorentz 
e 

term  -  c  x  B  of  course  vanishes.  This  approximation  is,  however, 

ma 


permissible  on  the  assumption  nD^  »  1  ,  and  the  "entropy  source  density" 
is  then  completely  due  to  the  so-called  collisions. 


For  the  time  variation  of  the  negative  entropy  within  a  closed  surface  £ 
we  get  the  expression 


dH 

so 

art 


/  •  (na  0  ln  fa)  dr  +  J  “TT”  dr 

v  *  v 

_  4  H 

/  na  c  ln  fa  dS  +  /  S*  dr  . 

£  v 


(2.15) 


If  the  first  contribution  disappears,  we  are  left  with  the  change  in 
due  to  collisions. 


dH. 

ao 


I* 


dr 


dcH« 

~ nr 


=  n. 


d  ln  f 
c  a 

dt 


(2.16) 


dH. 


Our  problem  is  now  to  show  that  collisions  can  only  lead  to 


“0  7  0. 


dt  < 

where  the  equality  sign  is  valid  when  the  state  of  equilibrium  is  attained. 


We  have  seen  earlier  that  collisions  in  ionized  gases  must  be  divided 
into  binary  or  short-range  collisions,  and  multiple  or  long-range  collisions. 
The  validity  of  the  H-theorem,  on  the  conditions  that  collisions  are  binary, 
was  first  proved  by  Boltzmann  [11]  .  More  unassailable  proofs  were  later 
given  by  Lorentz  [12]  ,  Jeans  [13]  ,  Fowler  [10]  and  others,  both  for  a 
simple  and  a  composite  gas. 

The  procedure  of  these  authors  is  substantially  that  described  by  Chap¬ 
man  and  Cowling  [4]  ,  and  consists  in  substituting  the  two-particle  colli  - 

d  f 

sion  integral  (1.26)  in  expression  (2.  16)  for  — ^  *  ,  and,  using  certain 
transformations,  showing  that  the  following  relation  is  valid, 


dcH« 

t r~ 


i  0 


(2.17) 


It  is  further  shown  that  the  equation  — -  =  0  and  detailed  balancing 

(i.e.  the  quantity  f#fp  is  unaltered  by  collisions  between  particles  of  the 
types  >  and  0  ),  can  only  be  valid  for  a  Maxwellian  distribution  function. 


In  these  treatments,  however,  the  effect  of  long-range  forces  of  ionized 
molecules  has  not  been  taken  into  account.  This  we  now  want  to  do. 


We  shall  show  that  the  H-theorem  remains  valid  provided  that  we  let 
the  molecular  interactions  be  represented  by  the  collision  term  on  Landau's 
form. 

Consider  an  ionized  gas,  composed  of  several  components  (  o  =  1,2, 
...)  with  distribution  functions  fa(r,0,t)  particle  masses  ma  and 
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d  f 
c  a 


particle  charge  ea  .  The  collision  term  ^  *B  the  sum  of  all  contribu- 
/  dcfa\ 

tions  I  — yj—j  p  ,  contributing  to  the  change  in  fa  from  collisions  with  par¬ 
ticles  0  .  According  to  ( 1 . 82) 


,  .  2  2 
2*Xe  e 
tt 


P 


m, 


"  .  (i  df'  f'  df  \ 

- f-  -  — S-)  dc'  ,  (2.18) 

do  \  mp  (jo  ma  / 


where  fjj  =  fp(o')  . 

The  total  negative  entropy  per  unit  volume  is 


2  2  v 

e«V  * 


H  =  =  /fa  In  f # d  c  . 

For  (  Vt°)  p we  get  *k«p =  2*  X 

=  -  !lM-/dclnfa-l./V.(-L 
\  6t  J  p  ma  a  do  \  mp  do' 

1  d  ln  f«  \ 

-  ST - ~  )  dc’  * 

ma  dc  y  * 


d  In  f' 


We  introduce  the  following  abbreviations, 

,  d  ln  f  ,  d  In  f 

%  =  %«=■  ■  sj-ir*-  •  sr— - 


and 


i.e. 


A.@  *  ®  v«p‘«fbdc 


(*) 


«P 


J  ln  £a  TT-  •  Aafl  dc  = 


P  '  ”5*7  '  77  '  *P 


k 

"m 


(2.  19) 


(2.20) 


/  4l  •  (In  faA.a)  dc  +  /  A  In  f  dc  . 

"a  do  *  m„  *  ap  4c  * 


(2.21) 


But  as  Agp  is  proportional  to  fa  (2.20),  the  first  integral  in  equation  (2.21) 
must  vanish  (assuming  that  fa  «♦  0  when  |e  |  »  ),  and  we  are  left  with 


(~rr)p  =  J  A«p (c) 


d  In  f 


do 


—  dc  , 


or 

m- 


(2.22) 


i  *  **  f. 


v  •  V«p(c,C). 
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For 


m 


we  get  in  a  completely  analogous  manner, 


(“Tt^a  =  kP«  Jdc/de-fpf.  •  V  •  Vpa(c,C). 


(2.23) 


The  value  of  the  integral  does  not  change  if  we  switch  about  the  primed  and 
the  non -primed  quantities, 

/  d  Hfl\  ,  din  f ' 

(-TJ1).  *  ‘f!.  STJ-  -  v  ■  V  («'•  •>  ■  <2  24> 

The  order  of  integration  may  also  be  changed, 

/  dH.\  ,  din  f' 

^-Wj«=  k0«J  doJdo' V«  i?6-1  V'C-c)  1 


3  k«fjJ  V  <«'•«>•  <2  “> 

By  summation  of  equations  (2.22)  and  (2.25)  we  get,  using  the  relation 
Vp#(c',c)  =  '▼ajJ<e»c')  .  (2-26) 


that 


(tt),  4  *  -k<*/dd  /  d«Vi>  T.P 


and  therefore  also, 
d  H  \  k 


(2.27) 

/  do  /  dc'  f#f^  Va#(c,  C)  •  V  •  V_(c,C)  .  (2.28) 


aa 


TJ  =  — L-  (u2  I  -  aa),  u  =  o  -  o' 


As 


the  expression  V  •  TJ  -V  has  the  form 


V  -  U-V  =  _L  [u2  ?-  I  •  V  -(?•■)(  a  •?)  ]  = 

u’ 


*  “T 

u 


[a2  ?2  -  (u  •  v)2]  5  o  ,  (2.29) 


and  the  expressions  (2.27)  and  (2.28)  are  always  negative.  They  equal  aero 
only  if  a  and  V  are  parallel  vectors . 
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We  have  thus  proved  the  H -theorem  even  in  the  following  case, 


d  H  d  Ha  d  H 

~TT  -  l  TT  ~  asp  "IT-  0  ’  *  «> 


*cH« 


+  dcHp 

TT  0  +  Tt  a 


<  0 


(2.30) 


These  results  must  be  interpreted  in  the  following  manner: 

The  entropy  for  each  gas  component  a  always  increases  as  a  result  of 
multiple  interactions  between  the  molecules  of  the  components,  and  attains 
its  maximum  value  when  Vaa  1 1  u  (2. 28), 


1 

m„ 


4  In  fj 
dc1" 


6  In  f 


IV  A  •  < 

- SL  =  J_u  =  —  (o  -  o') 

ma  do  K  1 


(2.31) 


for  all  e  and  o'  (  a  =  1,2,...) 
We  then  have 


1  dlnf« 


m. 


do 


1 


din  f 


«  .  1  . 


o  =  — - —  +  —  o 

ma  d  o1  *a 


and  consequently 

d  In  f 
_ c 

‘a  do 


1 

m. 


=  Const.  =  ~r—  C„„ 

*  ao 

<K 


<<=  -  °«o> 


(2.32) 


(2.33) 


with  the  solution 


m. 


=  Na  exP  t  '  177  (  0  "  °ao)  3 


eo' 


(«  =  1,2,  ...) 


(2.  34) 


The  entropy  for  each  gas  component  no  longer  grows  as  a  result  of  the 
interactions  between  the  molecules  of  the  gas  component  when  the  function 
of  state  has  become  Maxwellian  with  temperature  4a  and  mean  velocity 

°ao  * 

We  easily  verify  that  cao  and  are  really  the  mean  velocity  and 
mean  temperature  of  the  molecules  by  requiring  that  equation  (2.34)  satisfy 
the  conditions 


nB  =  /  f#dc  ,  n«e#0  =  Jvfgd*  * 

•/*»,(«■  °«o>2f«dc  •  <2-35) 

The  first  condition,  the  condition  of  normalization,  moreover  gives  the  fa¬ 
miliar  value  Nft  , 


Here  n  is  the  density  of  the  molecules  a  . 


(2.36) 
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We  further  see  that  the  entropy  will  no  longer  increase 
interactions  between  molecules  a  and  0  of  different  types 
V«P  11“’ 


as  a  result  of 
(2.27)  when 


i  ’lnfi 


nig  do' 


ma  cTc” 


-  i  -  -  i  <•  -  •’> 


(2.37) 


for  all  c 

and 

o'. 

We 

then 

have , 

1 

din  f 

a 

1 

0  = 

1 

d  In  f'p 

-  +  1 

m« 

dc 

mp 

do1 

* 

0 '  =  Const  =  —  C  , 

$  o 


(2.38) 


l.  e. 


I  d  in  f.  , 

mi -nr  -  j  ■ 


(2.39) 


which  has  the  solution 


m  . 

f«  =  N«exP('  zf  (°  ‘  *o>  ) 


(«  =  1.2.  •••)  • 


(2.40) 


Thus  the  equilibrium  distributions  are  once  more  Maxwellian,  but  in  this 
case  with  the  same  temperature,  <8  ,  and  mean  velocity,  CQ  ,  for  all 
components  of  the  gas  . 


Such  a  distribution  (2 . 40)  will  simultaneously  lead  to 
(for  all  a),  so  that 


0 


dH 

“3T  = 


0  . 


(2.41) 


Equation  (2.40)  will  thus  give  the  form  of  the  distribution  functions  in  the 
ultimate  state  of  equilibrium. 

In  a  gas  which  is  originally  far  removed  from  statistical  equilibrium  we 
must  consequently  expect  the  following  development: 

Encounters  between  identical  particles  will  initially  establish  local  equi¬ 
librium  distributions  for  each  gas  component,  each  with  its  own  temperature 
and  mean  velocity.  Next,  collisions  between  the  different  gases  will  lead  to 
equalization,  and  in  the  ultimate  state  of  equilibrium  the  components  will  have 
the  same  temperature  and  mean  velocity.  The  exact  time  development  can, 
of  course,  not  be  predicted  without  knowledge  of  the  characteristic  relaxa¬ 
tion  times  of  each  of  the  equalizing  processes  mentioned. 

These  relaxation  times  will  be  more  closely  studied  in  the  following 
section.  In  the  special  case  that  our  ionized  gas  is  a  plasma  consisting  of 
electrons  (ej.mi)  and  positive  ions  (eg.mg).  (mj  «  mg),  we  find  that 
the  proportionality  between  the  times  for  equalization  of  energy  between 

the  particles  a  and  {3  is 


"'ll  :  *22  :  *21  :  *12 


(2.42) 
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The  proportionality  between  the  equalization  times  for  momentum  is 


T 


1 

1  '* 


(2.43) 


This  makes  it  possible  to  consider  a  plasma  as  consisting  of  two  approxi¬ 
mately  independent  component  gases,  each  with  its  own  well-defined  tem¬ 
perature  equilibrium  is  established  for  each  component  much  quicker  than 
mutual  equalization.  This  is  especially  the  case  for  the  light  electrons. 


13.  CHARACTERISTIC  TIMES  FOR  THE  RELAXATION  PROCESSES 


In  the  preceding  section  we  showed  that  long-range  interactions  of  the 
type  that  will  occur  in  ionized  gases  will  make  the  system  approach  and 
attain  statistical  equilibrium. 


This  is  a  property  we  must  ascribe  to  all  possible  molecular  interac¬ 
tion  mechanisms  [10]  and  this  is  the  microscopic  explanation  of  the  macro¬ 
scopic  assumptions  that  the  entropy  must  always  increase  [Eq.  (3.84)]  . 


All  mechanisms  of  interaction  are  in  this  sense  equivalent  and  unimpor¬ 
tant  as  long  as  we  only  wish  to  describe  states  of  equilibrium.  If,  on  the 
other  hand,  we  wish  to  describe  the  development  of  states  that  are  not  in 
statistical  equilibrium,  the  types  of  interaction  are  highly  important,  as 
they  determine  the  relaxation  times  of  the  different  equalizing  processes. 

The  characteristic  times  Tap  and  for  equalization  of  energy  and 
momentum  between  two  groups  of  identical  molecules  a  and  0  may  in 


principle  be  evaluated  with  the  aid  of  the  known  expression  for 


P' 


The  average  changes  in  momentum  and  energy  due  to  encounters  for 
the  molecules  a  per  unit  time  is ,  according  to  the  transport  equation 
(1.57),  expressed  in  terms  of  the  velocity  Ca'  *  o#  “®ao 


Ra  =  ne  -jf<m«c'a)  >  Q«  *  na  -jf*  (im«  c'a2 )  •  (2-44) 

They  are  composed  of  the  contributions  from  collisions  with  the  single  com¬ 
ponent  gases, 


The  quantities  with  identical  indices,  a  =  p  ,  vanish  when  the  state  of 

(dc*a\ 

equilibrium  is  established  for  each  component,  i.e.  I  — j-g — -J  =  0  ,  while 

the  quantities  with  different  indices,  a  $  (3  ,  vanish  only  when  all  compo¬ 
nents  have  attained  the  same  temperature  and  mean  velocity. 


M 
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According  to  (1.57)  we  get  for  the  changes  in  the  mean  velocities  C( 
and  energies  because  of  collisions ,  ° 

ma  na  ST-  =  '  |  na  =  Qa  *  (2-46) 

If  we  can  now  find  expressions  for  and  Qft  it  is  possible  to  determine 
the  relaxation  times  for  equalization  of  mean  velocity  and  temperature  on 
the  basis  of  these  equations . 

Consider  a  gas  composed  of  electrons  and  positive  ions  with  masses 
mi  and  m2-  The  effect  of  the  collisions  may  be  represented  by  Landau's 
collision  integral  (1.82).  We  first  assume  that  each  component  gas  is  in  the 
state  of  equilibrium  at  different  values  of  cao  and  ,  so  that 


Ba  =  =  ^8  0  * 


«  ♦  0  • 


Approximate  expressions  for  these  quantities  may  be  found  by  taking 
into  account  that  the  mass  of  the  electron  is  much  smaller  than  the  mass  of 
the  ion,  mj  «  my  ,  and  that  the  thermal  velocity  of  the  electrons  is  much 
greater  than  that  of  the  ions,  provided  that  and  *2  are  °*  th*  same 

order  of  magnitude.  We  also  assume  that  the  difference  between  the  average 
velocities,  uQ  =  °io  ~  °20  ’  *B  8ma^  when  compared  to  the  electron 

velocity. 

The  tensor  U  in  the  collision  integral  depends  on  the  difference  between 
the  velocities  of  the  electrons  and  ions ,  and  may  then  be  developed  after  pow¬ 
ers  of  small  quantities.  In  the  case  that  the  ions  and  the  electrons  have  a 
Maxwellian  distribution,  we  get  according  to  Braginskij  [14]  , 


uc10  1 

(a)  =  B 12  «  uq  , 

^®20  1 

(b)  =  B2  =  B21  *  mjnj  uq  , 

Znl  “dT  ~  °1  =  °12  “  -Q21rR2l'  “o 


(2.47) 


(d) 


-  n. 


^2 

dt 


°2  =  Q2 1 


s  3 


mlnl 
m2  T1 


(  *1  "  *7) 


Here 


T1  * 


3lfm*, 


#/s 


T,  = 


4  YT7  n2 


(2.48) 


and  we  immediately  see  from  equation  (2.47a)  that  Tj  is  the  characteri¬ 
stic  time  of  change  of  momentum  of  the  electrons  because  of  collisions  with 
the  ions. 

The  characteristic  time  of  change  of  momentum  for  the  ions  because  of 
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collisions  with  electrons,  much  larger  (2.47b), 


m2n2 

mlnl 


1 


(2.49) 


Further,  equations  (c)  and  (d)  of  the  set  (2.47)  tell  us  that  the  relaxation 
times  for  energy  are 


12 


m, 
_  1  _ t 

-  *  », 


1  -  * 


m2n2 


21  1  mjnj  1  * 

i.e.  for  the  same  particle  density,  n^ 


=  n 


(2.50) 


2  ' 


we  have 


12 


=  T, 


1 

11  =  ^  T2  ~  $  mT 


1 


(2.51) 


The  characteristic  times 


Tir 


22 1 


1  2 

Tj  and  cannot  be  evaluated 


in  the  same,  relatively  simple  manner  because  a  series  development  of 

(  Va\ 

V  ~7t — /«  a^ter  small  quantities  is  no  longer  possible.  To  find  these  times 
we  must  evaluate  R  ma  and  Q  aa  in  other  ways  with  the  aid  of  equation 

(2.45).  Using  a  procedure  which  will  not  be  described  here,  Braginskij  has 
found  even  these  times, 

3 


% 


1  1 


11 


4  Ti  e  j  n  j 


m2  *2 


V. 


22 


and  moreover, 


4  V7  e  ^  n  ^ 


=  VT  T, 


■'fi 


(2.52) 


11 


=  T 


22 


(2.53) 


It  is  worth  noting  that  these  relaxation  times,  derived  with  the  aid  of 
Landau's  collision  integral,  agree  with  the  times  found  by  Spitzer  [15]  for 
the  establishment  of  equilibrium  of  energy  between  two  groups  of  particles 
with  approximate  Maxwellian  distribution, 


,  3m«md  /«.  .  ■»«  \  v» 

*e  ‘  TnTTgn-  lm«  ^e) 


(«,?  =  1.2).  (2.54) 
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and  with  Chandrasekhar's  [16]  times  tP  for  the  establishment  of  equili 
brium  of  momentum  between  two  groups  a  and  0  , 


tp  _  3  *2 


% 


i 


8*Ae2  e2 


Here  the  functions  4>(x)  and  G(x)  are  defined  by 

♦  <*>  =  -jb  /  ‘2  dt  ,  G(x)  =  ♦(-)  '  »*V)  . 

y«  ©  2x^ 


(2.55) 


(2.56) 


Spitzer's  times  (2.54)  for  “tga  agree  exactly  with  the  times  (2.50)  and 
(2.52)  if  we  put  mj  «  m^,  and  Chandrasekhar's  times  (2.56)  for  tP  agree 

with  the  times  (2.48),  (2.49)  and  (2.53)  on  the  same  conditions,  apart  from 
some  insignificant  numerical  factors. 


According  to  the  above  formulas  the  following  relations  exist  between 

the  relaxation  times  of  energy  and  momentum,  and  tP  ,  in  an  elec- 

,  _  \  op  o 

tron-ion  gas  (nj  s  n,)  , 


(2.57) 


and  moreover 
T11 


’!  * 


VT  Tj 


(2.59) 


The  characteristic  times  for  the  establishment  of  equilibrium  of  mo¬ 
mentum,  i.e.  isotropic  velocity  distribution,  thus  coincide  with  the  times 
for  the  establishment  of  equilibrium  of  energy  as  concerns  the  interactions 
electrons  ♦— *  electrons,  ions  «— »  ions,  and  electrons  ions.  But 

with  respect  to  interactions  ions  -*  electrons  we  see  that  the  ions  will 
make  the  distribution  of  the  electrons  isotropic  in  a  much  shorter  time  than 
they  will  need  to  equalize  the  temperatures.  We  must  thus  expect  the  elec¬ 
trons  to  acquire  an  isotropic  velocity  distribution  very  rapidly.  For  the 
ions,  on  the  other  hand,  the  two  processes  will  develop  with  equal  rapidity. 

It  might  be  interesting  to  make  a  comparison  with  the  results  obtained 
when  only  taking  the  short-range  collisions  into  account.  Chapman  and 
Cowling  give  the  following  value  for  the  collision  time  [4]  , 


T*3  =  «P 

ZfTxc1^  np** 


(2.59) 


where 


is  the  reduced  mass  for  two  particles  a  and 


P  * 


-  43  - 


mBm0 

Ha0  =  ,  while  0ffp  is  a  mean  impact  parameter. 

For  0ag  we  substitute  the  value  that  corresponds  to  a  relative  deflec¬ 
tion  if/ 2  , 


eaep 


rt  w“ 

°a  pwap 


(2.60) 


(wffp  is  the  relative  velocity  between  the  particles). 


m. 


m. 


If  we  take  into  account  that  Pjj  =  — —  ,  p^g 
(mj  «  mg) i  and  that  mean  values  of  wa^  are 


and  p12  z  mj. 


we  get  for  , 


(2.61) 
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(2.62) 
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The  mutual  proportions  are  thus  the  same  as  for  long-range  interactions 
(2.57). 


T 


1 

lc 


(2.63) 


But  we  immediately  see  that  the  deflection  time  for  short-range  collisions 
is  larger  by  a  factor  -  6  X  .  According  to  Section  4, 


k 


In  D/  q  «  In 


(2.64) 


The  quantity  k  grows  slowly  with  increasing  temperature  4  (=kT),  and 
decreasing  density  n  ,  and  stays  within  the  regions  of  temperature  and 
density  pertinent  to  plasma  physics  of  the  order  of  magnitude  -  20  [  15]  . 

The  collision  times  for  long-range  interactions  thus  are  less  than  the 
short-range  collision  times  by  a  factor  -  10^ .  A  grave  error  is  therefore 
committed  by  only  taking  short-range  collisions  into  account,  as  pointed  out 
by  Cohen  et  al.  [17]  .  This  is,  however,  permissible  in  the  case  of  gases 
with  low  temperature  and  high  density. 


11 
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14.  INTERACTIONS  WITH  THE  ENVIRONMENT 

In  Section  8  we  showed  that  the  force  term  <F_  F  >  of  the  averaged 

oam 

Liouville's  equation  (1.8)  may  be  split  into  the  following  parts, 

V"dc'+  J  i  r«e<V‘«Vdr'd»’  (2  65) 

The  first  part  is  the  external  macroscopic  force  field  while  the  second  part 
is  defined  as  the  internal  macroscopic  force  field, 

r«  *  |  /’,«?fidr'd<>'  •  <2-66> 

The  electromagnetic  field  involved  in  (2.66)  satisfies  Maxwell's  equations, 
as  shown  in  Section  3.  Finally,  the  last  part  of  (2.65)  is  due  to  the  corre¬ 
lations  between  the  different  kinds  of  particles  and  gives  origin  to  the  so- 
called  collision  term. 

The  important  thing  about  this  result  is  that  interactions  between  mole¬ 
cular  systems  in  the  correlation  region  (where  fag  -fa  f  =  f  „  4  0)  can¬ 
not  be  represented  by  macroscopic  force  fields.  "  " 

Part  of  the  macroscopic  field  Fa  is  due  to  the  external  field  1^0  and 
represents  the  effect  of  the  environment  on  the  system.  But  recalling  the 
above  statement,  we  cannot  expect  to  be  able  to  represent  the  total  effect  of 
the  environment  by  macroscopic  fields,  dependent  on  time  and  position. 

Apart  from  this  more  indirect  effect,  we  must  also  take  into  account 
the  possible  occurrence  of  direct  effects  having  the  character  of  random 
collisions  between  gas  molecules  and  the  molecules  of  the  environment,  the 
mutual  distance  of  which  is  less  than  a  certain  correlation  length.  The  en¬ 
vironment  may  consist  of  other  gases ,  fluids  or  solid  matter .  Thus ,  in 
general,  a  knowledge  of  the  detailed  microscopic  structure  both  of  the  system 
and  the  surroundings  is  absolutely  necessary  to  give  an  exact  description  of 
the  development  of  the  system  by  equation  (1.62).  This  really  means  that  we 
extend  the  system  to  include  part  of  the  surroundings. 

This,  however,  we  do  not  wish  to  do,  and  we  must  consequently  here¬ 
after  neglect  direct  "collisions"  between  system  and  environment  and  assume 
that  the  total  effect  of  the  surroundings  on  the  system  may  be  represented  by 
external  macroscopic  force  fields  which  depend  on  time  and  position: 

F  =  F  (r,t).  This  assumption  is  justified  if  the  gas  is  sufficiently  far 

B  O  8  0 

away  from  the  environments  or  if  the  molecules  at  the  surface  of  the  gas  are 
in  local  statistical  equilibrium  with  the  molecules  in  the  nearest  environment. 

If  the  environment  is  constituted  by  another  (ionised)  gas ,  we  have  seen 
that  eventually  the  collisions  will  invariably  lead  to  the  attainment  of  the 
equilibrium  state.  We  must  expect  similar  conditions  when  the  surroundings 
consist  of  solid  matter  or  fluids .  (The  entropy  of  the  system  can  admittedly 
diminish  because  of  interaction  with  the  environment.  This,  however,  is  al¬ 
ways  accompanied  by  a  corresponding  increase  in  the  entropy  of  the  environ¬ 
ment,  so  that  the  total  change  is  positive.) 
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15.  EXTERNAL  TIME -DEPENDENT  MAGNETIC  FIELD 


Before  leaving  the  study  of  the  statistical  properties  of  ionized  gases, 
we  shall  consider  a  state  of  affairs  which  at  first  appears  paradoxical,  but 
which  actually  has  a  natural  explanation. 

We  have  seen  that  the  only  physical  effect  which  contributes  to  the 
change  in  entropy  when  the  "entropy  flow"  (na  c  In  fa  )  is  zero,  is  the  so- 

d  H 

called  "entropy  source  density"  — .  The  latter  is  always  negative. 


As  a  result  of  molecular  interactions  the  entropy  of  a  system  will  de¬ 
crease  towards  a  lower,  stable  value  (Ha  )  while  the  state  of  the  system, 

represented  by  fa,  will  simultaneously  approach  statistical  equilibrium , 

i.e.  iaotropto  Maxwell -Boltzmann  distribution. 


Suppose  that  we  have  an  external  magnetic  field  in  the  gas,  and  that 
there  is  present  a  time  change  of  the  field  with  such  a  short  period  x  that 
the  collisions  are  unable  to  "smooth  out"  the  effect  in  less  then  the  time 

2  %  eB 

t  (t  «  t  ).  On  the  condition  -  «  t  (w  -*  - )  the  magnetic  moment 

r  r  w  me  " 

1  2 
t  me 

of  a  particle,  p  =  - — = —  ,  is  an  adiabatic  constant,  and  the  transverse 

B 

particle  energy  will  therefore  change  proportionally  to  B  ,  while  the  longi¬ 
tudinal  one,  |  mc^  j  ,  is  of  course  unchanged..  Analogous  conditions  exist 

even  when  ?-5-  «  X.  is  not  satisfied, 
w  K 


It  appears  as  if  the  system  retreats  from  statistical  equilibrium  because 
of  the  introduction  of  anisotropy.  If,  however,  we  have  a  look  at  the  expres- 
dH 

sion  for  — ,  we  find  that  there  is  no  corresponding  change  in  the  entropy. 


as  the  term  na  ■  Fa  In  f  „  (2.  15), 

dH_ 


from  external  fields,  disappears. 


representing  the  contribution  to 
The  solution  of  this  apparent  para¬ 


dox  is  immediately  seen  if  we  let  the  magnetic  field  revert  to  its  original 
value.  Provided  that  the  collisions  have  not  had  the  time  to  cause  disturb¬ 
ance,  we  shall  again  have  the  original  isotropic  state.  The  process  is  re¬ 
versible,  at  least  as  far  as  the  system  gas -magnetic  field  is  concerned, 


and  is  consequently  not  accompanied  by  a  change  in  entropy 


A  periodic  change  in  the  field  with  characteristic  time  equal  to  or 
greater  than  the  relaxation  time  for  the  gas  (x  Z  xr),  however,  will  result 
in  decreasing  entropy,  as  the  collisions  then  have  the  time  to  smooth  out 
of  the  difference  between  the  transverse  and  the  longitudinal  energy  [18]  . 

Ultimately  then,  only  the  collisions  can  lead  to  a  change  in  entropy, 
but  any  (reversible)  deviation  from  the  isotropic  Maxwell  distribution  in 
the  state  of  the  system  will  enable  the  collisions  to  increase  the  entropy 
of  the  system  further. 
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Finally  a  word  about  the  concept  of  "collisions".  The  "collision  term" 
in  Boltzmann's  equation  has  been  the  source  of  much  confusion,  evidenced, 
for  instance,  by  the  many  forms  of  the  term.  This  is  the  more  strange 
when  one  considers  the  fundamental  importance  of  the  so-called  collisions 
for  the  statistical  development  of  the  system. 

This  untenable  situation  is  clarified,  however,  when  we  take  our  point 
of  departure  in  the  exact  Liouville  equation  and  not  in  the  Boltzmann  equation 
with  a  more  or  less  arbitrary  collision  term. 


16.  THE  EQUILIBRIUM  CONFIGURATION  FOR  THE  SYSTEM 

16.1  DEDUCTION  FROM  MICROSCOPIC  THEORY 

Statistical  equilibrium  for  a  system  is  characterized  by  the  fulfillment 
of  two  requirements: 

(1)  Internal  equilibrium  for  the  system. 

(2)  Equilibrium  between  the  system  and  the  environment. 

A  macroscopic  system  always  has  an  internal  structure,  and  can  conse¬ 
quently  be  subdivided  into  component  systems,  i.e.  molecules,  groups  of 
molecules  etc. 

The  requirements  (1)  and  (2)  are  again  valid  for  all  macroscopic  compo¬ 
nent  systems  in  statistical  equilibrium.  In  this  case,  however,  the  environ¬ 
ment  with  respect  to  a  single  given  component  system  includes  all  other  com¬ 
ponent  systems . 

We  must  assume  that  all  (component)  systems  have  a  natural  tendency 
toward  statistical  equilibrium  because  of  the  intermolecular  interactions  in 
the  system  and  with  the  environment.  If  we  denote  the  negative  entropy  of 
the  (component)  system  by  H  ,  we  can  express  the  irreversible  tendency 
in  the  following  manner, 


sytt. 


•nv. 


(2.70) 


Both  interactions  within  the  system  and  with  the  surroundings  will  lead  to 
increasing  entropy.  The  requirements  (1)  and  (2)  can  now  be  expressed  as 


(1) 

(2) 


=  0  , 


(2.71) 


The  relations  (2.70-71)  are  only  valid  for  maoroaoopto  component 
systems.  For  microscopic  or  molecular  systems  the  so-called  "principle 
of  microscopic  reversibility"  is  applicable,  meaning  that  from  a  micro¬ 
scopic  point  of  view  the  development  is  invariant  with  respect  to  inversion 
of  time  (t  -*  -  t).  (In  the  case  that  we  have  an  external  magnetic  field,  its 
direction  must  be  reversed  (  B  -*  -  B  ). )  The  equations  of  motion  of  the 
single  molecules  are  symmetrical  with  respect  to  time. 

In  the  preceding  we  have  treated  a  composite  gas  partly  as  a  system 
and  partly  as  being  composed  of  independent  component  systems  (gases). 
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In  the  former  case  we  defined  the  thermal  velocities  in  relation  to  the  mass 
velocity  c0  for  the  whole  gas  (  Ca  =  o0  -  0o),  while  in  the  latter  case 
the  thermal  velocities  are  more  naturally  defined  in  relation  to  the  mean 
velocities  caQ  (Ca’  =  °a  "ca0)  of  the  component  gases.  In  this  case  we 

further  define  average  values  and  flux  vectors  of  molecular  quantities 
for  each  component  gas  in  the  same  way  as  done  earlier  for  a  simple  gas. 
We  then  get  as  many  sets  of  such  quantities  as  there  are  components  in  the 
gas,  and  as  many  sets  of  macroscopic  transport  equations. 

Under  the  first  point  of  view  the  interactions  between  the  component 
gases  will  have  the  form  of  internal  forces,  the  effects  of  which  cancel  when 
the  gas  is  viewed  as  a  whole.  Under  the  latter  point  of  view  the  interactions 
between  one  component  gas  and  the  other  components  will  have  the  form  of 
forces  between  one  component  and  the  surroundings.  The  requirements  (1) 
and  (2)  will  be  satisfied  in  both  cases  if  the  conditions 


P  ‘ 


0  , 


(2.72) 


for  all  a  and  0  ,  are  satisfied.  The  first  condition  means  that  the  distri¬ 
butions  fa  have  the  form  (2.40), 

e  ,  o=l,2 . 

(2.73) 


The  mean  velocities  0aQ  and  the  temperatures  4a  are  the  same  for  all 
component  gases,  i.e.  0Q  and  £  .  The  second  condition  is  satisfied  if 

we  disregard  possible  collision  terms  representing  interaction  with  the  en¬ 
vironment  in  Boltzmann's  equation.  When  we  impose  the  stationary  condi¬ 
tion  ( -gy.  =  0  )  this  equation  gets  the  form 


df 
_ a 

dr 


+  +  ca  x  wa) 


H* 

do. 


=  0  , 


(2.74) 


or,  if  we  introduce  the  thermal  velocity  C„  =  c  -  0^  and  divide  by  f  , 

ft  ft  O  7  ft 


d  In  f  d  In  f  d  In  f 

-  +  C  .  -  +  [F  -  c  •  JL  c  +  c  x  h  ]  ■  - 5 —  + 

dr  0  dr  1  a  °  dr  °  o  oJ 


din  f 

+  (cax«j. 


d  In  f 


dC„  dC„ 

a  a 


—  C#  :  -1-  co  =  0  . 

*  dr  o 


(2.75) 


In  contrast  to  earlier  conventien  Fa  does  now  not  include  the  magne- 

e  aB 

tic  force  c  x  o»  ,  w_  =  =r~—  • 
a  a  #  m  c 


From  this  equation  we  can  now  derive  the  equilibrium  conditions  for  the 

o'  '  ’  • a  ~a 


macroscopic  quantities  na ,  oQ,  4,  Fa  and  . 


It 
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d  In  ffl 

In  equation  (2.75)  we  substitute  the  expressions  obtained  for  - 

din  f  dr 

and  - - -  from  equation  (2 . 73),  getting 
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+  C«  Ca 
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dr 


mg  d  In  ♦  s 
2  *  dr 

(2.76) 


This  is  an  identity  in  Ca  ,  and  the  coefficients  of  the  powers  of  the  third, 
second,  first  and  zero  degree  must  therefore  be  zero.  This  gives  the  equa¬ 
tions  (2.77),  (2.78),  (2.80)  and  (2.85). 

We  get,  to  start  with, 
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i.  e. 


d  In  > 

dr 


0  i.e.  $  =  Const.  (2.77) 


Further,  using  (2.77), 
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(2.78) 


The  deformation  tensor  e  is  thus  zero,  and  the  gas  consequently  moves 
in  a  screw  motion  like  a  rigid  body. 


o  =  u  +  m  x  r  , 
o 


(2.79) 


where  u  and  *  are  constants. 

We  further  derive  from  equation  (2.76) 


'  d  In  n  m_ 
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As  o0  has  the  form  (2.79),  the  term  oQ  •  cQ  may  be  transformed 

into  m  x  cQ  suid  equation  (2.80)  can  be  written 


d  In  n  IB 
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a.  ■  •>„  +  « 


Now, 


c  x»  =  («  x  r  )  x  «  s  w  2  r  -  ( •*  •  r  )  •*  =  w2  r  x 


(2.81) 


(2.82) 
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where  lj_  ie  the  component  of  r  normal  to  the  axis  of  the  screw  motion 
(2.79).  We  can  therefore  write  oQ  x  as  as  the  gradient  of  a  centrifugal 
potential  $c> 

eox-=  4>c  =  -i<**Zr2  •  (2.83) 


Substitution  in  equation  (2.81)  shows  that  it  must  be  possible  to  derive  the 
vector  Fa  +  CQ  x  ug  from  the  potential  pg  , 

F«  +  co  x  *«  =  ‘  •  4>a  =  '  ln  n«  -  <j>c  •  <2-84) 

Therefore, 

■ 

n  =  n  e  ,  n  =  const.  (2.85) 

8  a  o  80  '  1 


The  last  relation  obtained  from  equation  (2.76), 


d  ln  n 


0 


(2.86) 


means  that  the  force  Pa  must  everywhere  be  normal  to  the  velocity  0Q  ,  so 
that  no  work  is  done  and  the  field  is  "force-free".  The  velocity  cQ  must 
lie  along  an  equi -potential  surface  =  Const. 

In  the  case  of  a  pure  translational  motion  0Q  *  a  ,  the  potential  is 
thus  cylindrical,  pa  =  ^(r,>)  ,  with.the  cylinder  axis  (z-axis)  along  «, 
while  for  a  rotational  motion  e0  =  is  x  r  the  potential  is  rotationally  sym¬ 
metrical  with  it  as  the  axis  of  rotation:  *  $a(r,z).  If  both  ■  and 

•  ♦  0  ,  then  p  is  constant  along  helical  trajectories  with  O  ( ||  at )  as  the 
axis. 

«  ea 

The  contribution  to  ra  from  an  electric  field  E  is  *_  =  - B 

a  mB 

If  6a  denotes  the  rest  of  Fg  (gravitational  fields,  etc.)  we  have, 

Fa  =  *a  +  °u  (2.87) 


Provided  that  Ga  may  be  derived  from  a  potential 


d$ 


°«  - 


8G 


dr 


equation  (2. 84)  gives 


*8  +  co  *  **i»  =  -~(4>.  ‘  $aG)  . 


dr 


(2.88) 


(2.89) 


meaning  that  the  vector  is  vorticity-free,  or,  as 
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t 


a 


and  m 

a 


l 

c 


we  have 

1i-x(E+iooxB)  =  0.  (2.90) 

But  this  condition  means  that  the  magnetic  flux  through  a  surface  that  follows 
the  motion  of  the  gas  will  be  constant  during  the  motion  [15]  ,  and  the  mag¬ 
netic  field  will  consequently  be  "frozen  into"  the  gas  when  the  latter  is  in 
statistical  equilibrium. 

According  to  the  above,  the  assumed  tendency  of  the  field  to  remain 
frozen  in  [15]  is  a  consequence  of  the  natural  tendency  of  the  system  to¬ 
wards  equilibrium  (Section  12). 

The  equations  (2.77),  (2.79),  (2.85)  and  (2.86)  as  well  as  the  two  addi¬ 
tional  conditions  that  the  temperature  and  the  mean  velocity  be  the  same  for 
all  component  gases  a  ,  express  the  characteristic  conditions  of  a  gas  in 
statistical  equilibrium, 


i  =  Const.,  c=o+»xr,  n  *  n  e 


-  -r-  <  V  ) 


d<k 


=  0  ,  , 


ao 


O  =0 
ao  o 


(2.91) 


°  d 

The  three  first  equations  may  be  expressed  on  differential  form 


din  A 
dr 


0  , 


e  =  o  , 


dr 


=  0  , 


(2.92) 


d  M 

when  - —  is  defined  by 

dr 


dn„  In  n  m 

=  - 5 - —  (F  +  c  xa  ),  a 

dr  dr  *  0  °  «  o 


+  m.  (2.93) 


In  an  arbitrary  state  these  three  quantities  constitute  a  measure  of  the  devia 
tion  from  statistical  equilibrium.  In  that  case  «Q  x  *»  must  be  replaced  by 


It  follows  from  the  set  (2.91)  that  J  *  =  0  (equation  (1.47)),  and  the 
current  density  is  Mien  equal  to  the  convective  flow  Qe  oQ  .  Maxwell's 
equations  therefore  assume  the  form  (^-  =  0), 


4  *  Q 


-k  *  -  0 


(2.94) 
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This  set  of  equations,  together  with  equations  (2.80),  (2.85)  and  (2.79), 


E+  J  on  xB  =  -  ^  , 

co  dr 


m 

(♦■  —  ♦«)  . 

ea 


(2.95) 


c  =  0  ,  oQ  =  u  +  »  xr,  «e  =  2  eana  * 

a 


gives  the  field  configurations  that  are  characteristic  of  statistical  equili- 
brium.  We  must  assume  that  the  fields  possess  a  natural  tendency  towards 
this  equilibrium  configuration. 

We  have  here  assumed  that  the  only  forces  are  electric  or  magnetic 
ones  (Ga  =  0 ). 

,  The  solution  of  the  equations  (2.94)  are,  when  disregarding  border 
effects , 


E=._1_A+E.  B=-i- 

dr  *  o  d  i* 


4TXA*  “o 


91  d  t 


♦.<*>  'l-fal  ■  k-  4-fr^ 


(2.96) 


Here  B  and  Bq  are  the  external  fields  that  satisfy  the  homogeneous 
Maxwelr  equations . 


In  the  special  case  that  CQ  =  Const.  =  a  we  derive  the  following 
relation  for  the  internal  fields  (exclusive  of  E0  and  B0)  from  equations 
(2.96), 


B  .  I  E  X  c  . 

C  0 


(2.97) 


Non-relativistically ,  the  internal  magnetic  field  must  therefore  be  zero, 
and  in  this  case  the  external  magnetic  field  is  frozen  in  the  gas . 

In  the  case  of  a  rotating  gas  the  conditions  become  more  complex.  The 
relation  (2.97),  for  instance,  is  no  longer  valid. 

The  above  equilibrium  conditions  (2.91)  and  (2.92)  are  microscopically 
founded  and  rest  on  the  validity  of  the  H -theorem.  In  classical  thermodyna¬ 
mics  the  equilibrium  conditions  are  derived  by  macroscopic  arguments  that 
have  the  character  of  postulates.  It  might  be  interesting  to  see  how  the 
usual  postulates  must  be  modified  to  be  consistent  with  the  results  of  micro¬ 
scopic  theory. 


16.2  COMPARISON  WITH  MACROSCOPIC  THEORY 

The  thermodynamical  expression  for  the  entropy  (per  unit  volume)  for 
an  ideal  gas  (nD^  »  1)  is  [4]  . 

Sa  =  -  nBk(ln  na  +  yin  ^  gfcy  "7^  •  (2.98) 

while  the  enthalpy  Wa  is  given  by 


It 
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(2.99) 


W. 


T  n«  kT 
2  * 


The  thermodynamic  or  Gibbs  potential  Ga  =  Wa  -ST  is  consequently 
given  by  ( &  =  kT). 

,  m 

Ga  =  na  »  (In  na  +  i  In  +  1)  .  (2.  100) 

Further,  the  chemical  potential  p°  is  defined  as  thermodynamic  potential 
per  molecule,  * 

,  m 

I*®  »  G«/na  •  *(lnn#  +|  In  jJL  +  1)  .  (2.101) 

In  the  case  that  we  have  external  and  internal  fields  of  force  Fa  that  may  be 
derived  from  a  potential  ma  the  expression  (2.  101)  must  be  replaced  by 

+  ma4a  .  (2.102) 


The  conditions  that  are  usually  assumed  for  thermodynamic  equilibrium 
are  constant  temperature,  velocity  and  chemical  potential, 

£  =  Const.  ,  oQ  =  U  ,  |»  =  Const.  (2.  103) 

These  three  conditions  may  be  written  on  differential  form. 


din  & 
dr 


=  O  , 


where 


dr 


din  n. 
a 

dr 


maF 


o  * 


dr 


0. 


(2. 104) 


(2. 105) 


These  are  the  conditions  assumed  by  for  instance  H.  S.  Green  [8]  .  He 
further  assumes  arbitrarily  that  the  relation 

B  =  J  E  x  c0  (2.  106) 

is  valid  in  equilibrium  between  the  internal  fields.  Non-relativistically  we 
then  have  B  =  0  . 

The  conditions  (2. 103)  or  (2.  104)  and  (2.  106),  which  must  be  considered 
as  arbitrary  macroscopic  postulates,  in  the  case  of  no  rotation  (  «s  =  0  ) 
and  no  magnetic  field  (  «u  =  0  )  ,  then  are  quite  equivalent  with  the  three 
first  conditions  (2.91)  or  (2.92)  which  are  founded  on  microscopic  theory. 

By  the  latter  procedure  we  get  in  addition  the  three  last  of  the  conditions 
(2.91),  which  cannot  be  macroscopically  proved. 

Consider  the  case  that  the  body  has  a  constant  rotation  •*  .  If  we  then 

2  2 

add  the  centrifugal  potential  m  a  |  mfl  m  r  to  the  chemical  po¬ 

tential  (2.102), 
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Ha  =  +  ma$a  -  i  ma  w2r2  ,  (2.107) 

and  further  require  that  the  magnetic  force  eQ  x  m  cam  be  derived  from  a 
potential  which  may  be  included  in  $a,  i.e. 


d 

dr 


x  (c  x  B  ) 

'  o  ' 


0 


(2. 108) 


(B  is  the  sum  of  the  internal  and  external  fields),  we  obtain  equilibrium 
conditions  that  are  quite  equivalent  with  the  conditions  (2.91)  or  (2.92). 
The  three  last  conditions  (2.91)  must,  however,  be  separately  postulated, 


o  • 
o 


0  , 


c 

to 


c 

o 


(2.109) 


We  have  seen  that  the  conditions  (2.91)  and  (2.92)  are  sufficient  to  se- 
cure  that  the  entropy  production  is  zero  and  that  the  state  is  stationary.  We 
shall  therefore  use  these  conditions  as  definitions  of  thermodynamic  or  sta¬ 
tistical  equilibrium.  We  must,  however,  be  aware  of  the  apparent  absence 
of  a  general  agreement  as  to  which  equilibrium  conditions  are  the  correct 
ones . 

The  first  part  of  this  section  is  essentially  a  generalization  of  the  results 
of  Chapman  and  Cowling  for  simple  gases  without  electromagnetic  fields  [4] 
to  the  electromagnetic  case  for  composite  gases. 
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PART  III.  PHENOMENOLOGICAL  THEORY  OF 
TRANSPORT  PHENOMENA. 

17.  TRANSPORT  PHENOMENA 

The  denotation  "transport  phenomenon"  is  applied  to  certain  phenomena 
occurring  in  a  gas  which  is  not  in  statistical  equilibrium.  We  have  seen 
that  any  deviation  from  equilibrium  will  set  off  microscopic  equalization 
processes,  which,  in  accordance  with  Le  Chatelier's  principle,  seek  to 
counteract  their  cause  and  lead  the  system  back  to  the  equilibrium  state. 
The  macroscopic  counterpart  of  the  microscopic  equalization  tendencies 
are  the  so-called  transport  phenomena.  The  fact  that  a  system  is  not  in 
statistical  equilibrium  can  always  be  traced  back  to  disturbing  interaction 
with  the  environment,  with  the  result  that  the  equilibrium  conditions  (2.91) 
can  no  longer  be  satisfied.  (We  then  disregard  possible  changes  in  the  in¬ 
ternal  structure  of  the  molecules.)  These  disturbances  set  up  gradients 
(possibly  time  variant)  in  temperature,  velocity  and  density,  that  do  not 
agree  with  (2.92).  For  composite  gases  the  equilibrium  between  the  com¬ 
ponent  gases  will  moreover  be  destroyed.  We  thus  get  phenomena  like  dif¬ 
fusion,  electric  current,  transport  of  momentum  and  heat  conduction. 

The  microscopic  explanation  of  this  is  that  the  molecules ,  during  their 
disorganized  thermal  motion,  will  "transport"  their  molecular  properties 
(e.g.  number,  electric  charge,  momentum  and  energy)  from  place  to  place. 
In  statistical  equilibrium  this  transport  is  completely  random,  so  that  the 
microscopic  mean  values  disappear.  Deviation  from  equilibrium,  on  the 
other  hand,  means  that  this  is  no  longer  the  case. 

Let  $  =  p  (r,  C  ,  t)  denote  the  molecular  property  as  before.  Then 

n$  =  I  n#  =  I  /♦afadCa  (3.1) 

a  a 

is  the  corresponding  macroscopic  property.  The  transport  of  this  property 
is  given  by  the  so-called  flux  vector 

n?C  =  I  n0  f0Ca  =  Z  J  ♦aCafadC0  ,  (3.2) 

which  gives  the  flux  of  the  property  $  through  a  surface  which  follows  the 
mean  motion  of  the  gas . 


17.1  DIFFUSION  OR  NUMBER  TRANSPORT 

If  we  assign  to  $  the  molecular  property  "number",  i.e.  ^  =  1, 
we  get  from  (3.  1), 

n«C«  -  /CafadC«  •  (3-3) 

This  vector  gives  the  number  flux  of  molecules  a  in  relation  to  the  mean 
motion  of  the  gas.  Another  name  for  (j a  tbe  "diffusion  velocity"  of 
the  molecules  a  .  If  we  have  two  component  gases  (  a  =  1  and  2),  we 
are  usually  interested  in  the  relative  diffusion  velocity. 


■|-*2  ■  ST 
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JC2f2dC2  = 

=  ®i  "  ®2  =  °10  ‘  °20  =  Uo  ’  (3-4) 

which  equals  the  difference  between  the  mean  motions  of  the  two  gases. 

We  have  a  special  case  when  the  molecules  of  the  gases  1  and  2  are 
identical.  The  vector  C  ^  '  C  j  then  gives  the  "self-diffusion"  of  the 

molecules,  i.e.  the  relative  diffusion  of  two  different  groups  of  identical 
molecules . 

The  diffusion  is  evidently  zero  when  the  gas  is  in  statistical  equilibrium 
(Eq.  (2.73)),  naCa^°^  =  0  .  The  symbol  ^  denotes  statistical  equili¬ 

brium  . 


17.2  ELECTRIC  CHARGE  TRANSPORT 
The  total  electric  current  density  is  given  by 

i  =  J  eftca  =  (  2  naea)  °Q  +  2  a  .  (3.5) 

and  is  the  sum  of  the  "convection  current"  ((  J  n  e_  )  c  )  and  the  other  term 

■  o  o 

which  gives  the  proper  conduction  current  J  *  .  The  convection  current  is 
zero  when  the  gas  is  electrically  neutral. 

The  flux  vector  for  electric  charge  is  obtained  from  equation  (3.2)  by 
putting  $  =  e  , 


neC  =  2  nae#CB  =  J  J  e«  {adCa 


(3.6) 


This  is  identical  to  the  conduction  current  j  *  , 

J*  =  neC  =2  naeaCa  =  2  J*  #  . 


(3.7) 


The  contribution  to  this  from  the  molecules  a,  i.e.  J  *  =  nae05ft» 

thus  is  proportional  to  the  diffusion  vector  n#Ca  • 

For  a  binary  gaB  mixture  we  get  from  equation  (3.7)  and  the  relation 


Jnam«Ca  =  0  , 


n.n,(m,e.  -m.e,) 

j*  ■  12  V  12  <e,-e2)  . 


(3.8) 


In  this  case  the  electric  flux  vector  is  proportional  to  the  relative  diffusion 
velocity.  If  in  a  plasma  we  denote  the  electrons  by  index  1  and  the  ions  by 
index  2,  i.e.  m^/m^  «  1,  equation  (3 . 8)  will  give  approximative^, 


J*  =  nje j(  g  j  -  C2) 


(3.9) 


14 
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i 


The  current  density  is  zero  in  statistical  equilibrium ,  (Eq.  (2.73)), 

•<«>  =  0  . 


17.3  TRANSPORT  OF  MOMENTUM 

We  have  previously  shown  (Section  6)  that  the  flux  tensor  for  momentum, 
or  the  stress  tensor,  is  found  by  putting  ■  mC  in  equation  (3.2).  We  then 
get, 


p  =  nmCC  =  I  n8maCaCa  =  I  pa  (3.10) 

The  expression  can  be  transformed  by  introducing  the  total  velocity 


p  =  I  namaCa  Ca  +  Q  CQ  Cq  .  (3.11) 

In  statistical  equilibrium,  where  the  distribution  is  given  by  (2.73),  we 
easily  find  for  p  , 

p(o)  =  pi,  P  =  I  na  =  n*  =  nkT  ,  (3.12) 

a 

as  $a  =  ♦  .  All  shear  forces  disappear,  and  the  stress  tensor  repre¬ 

sents  a  purely  hydrostatic  pressure.  According  to  equation  (1.54)  this  can 
also  be  expressed  by  the  condition 

J}(°)  _  p(o)  .  p  i  =  ©  .  (3.  13) 

17.4  HEAT  TRANSPORT 

We  have  earlier  found  that  the  energy  flux  vector  is  obtained  by  substi¬ 
tution  of  $  =  E  =  %  mC^  in  equation  (3.2).  The  result  is 

q  =  nimCZC  =Jna*m«c2«C«  =**«  •  (3-  14) 

The  heat  flow,  as  well  as  the  diffusion,  conduction  current  and  shear  forces, 
disappear  when  fa  is  isotropic  in  the  velocity  Ca  .  This  is  especially  the 

case  for  Maxwell -Boltzmann  distribution  (2.73)  in  statistical  equilibrium, 
i.e.  q  (°)  =  0  . 


18.  THE  TRANSPORT  EQUATIONS 

In  Section  7  we  derived  the  macroscopic  transport  equation  for  the  quan¬ 
tity  $  where  $  is  a  molecular  property  (equation  (1 . 57)).  Three  molecular 
properties  are  of  special  interest.  They  are  the  so-called  summational  in¬ 
variants 

$  =  i  mC2  , 


4>  =  1  . 


4>  =  m  C  , 


(3.15) 
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that  is  to  say  molecular  number,  momentum  and  energy.  The  reason  that 
these  three  quantities  are  in  a  class  by  themselves  is  that  the  total  amount 
of  these  properties  is  conserved  during  collisions  between  molecules, 


dc$ 

~TT 


=  I 
a 


«c*. 

71 


=  I 
a 


dcf« 

jr 


dC. 


=  0 


(3.16) 


One  may  satisfy  oneself  that  equation  (3.  16)  is  correct,  when  — — —  is 

d 1 

Landau's  collision  term,  by  a  direct  computation  which  is  given  in  the 
Appendix.  Equation  (3.  16)  is  also  valid  for  the  familiar  two-particle  col¬ 
lision  integral (1.26). 

The  quantities  (3.  15)  are  the  only  ones  having  this  property. 

For  a  simple  gas  the  transport  equations  take  the  form 


(a) 

dn 

ar  +  n 

o 

II 

o 

o 

9 

(b) 

dCo 

«-ar- 

-  QP  +  — 
w  dr 

•  JP  =  o  , 

(3.17) 

dU  + 
~3T  + 

U  -4-  •  c  + 
dr  ° 

■JT'  "  ■,e 

■  F  *  *  ;  -fr  °o  ■  °- 

2  3 

where  U  =  nJmC  =  £  n$  is  the  average  molecular  energy,  and  ]p  and 

q  are  the  stress  tensor  (3.  10)  and  heat  flow  (3.  14),  respectively.  These 
three  equations  thus  express  the  conservation  of  molecular  number,  momen¬ 
tum  and  energy  in  the  gas . 

For  a  gas  with  several  components  ( a  =  1,2,  .  .  .)  we  obtain  similar 
sets  for  the  components , 


(a) 

(b) 


dn  A  .  _  4  1 

-ar+n«  « '  °o +  ' (n*  °«)  *  “«  tt  =  0  ’ 


-TT-*".®.  eo*-k-  -«.<* 


d  C 


e'ra  -~ar>+  9«  V  Tr* ° 


=  n . 


dm  C» 
c  a  a 


dU  de 

"ar~+  u«  ir  '  °° +  dr'  q« "  9«c«‘  (F« '  "af)  + 


dr 

+  :  -k  °o  =  n, 


dc  i  m«C 


dt 


(3.18) 


By  adding  the  corresponding  equations  for  all  component  gates  and  using 

equation  (3.  16)  and  the  relation  Z  Q  C  =  0  we  derive. 

.88 
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(a) 

dQ  , 

ar  + 

d 

a?  • 

O 

o 

II 

o 

(b) 

0dco 

Q-ar 

dU  + 
ar  + 

•  P  -  1  SffF  a  =  0  , 

a 

(3.19) 

(c) 

ui. 

dr 

c  +  -i-  •  q  -  2  •  F 

°  dr  a  a  a  a 

+  J>  :  -2- 

“  dr 

0o  3  0  • 

Here 

U  =  2  ua  = 

j  n-»  ,  p  =  2  pa  and 

“  ft 

are  the  average 

molecular  energy  per  unit  volume,  the  atreaa  tensor  and  the  heat  flow  for  the 
whole  gas,  respectively. 

In  the  case  that  only  electric  and  magnetic  force  fields  are  present, 
e. 

F. 


m 


01  (E  +  jj-  ca  x  B)  ,  the  equation  (3.  19)  can  be  further  transformed. 


(a) 

(b) 

(c) 


+  Q  •  C  =  0  , 
dT  M  <jr  o 

d  c 

o 

Q  -ar- 


-  Qj*+  ±  coxB)  J*x  B+  P  =  0  ,  (3.20) 


+  ti.q.J*.(B+IcnxB)+p:l  c, 

at  dr  °  dr  M  c  o  dr  { 


=  0  . 


Here  Qfi  =  I  nfl  ea  *  I«e#  and  J*  =  2  na  e  a  Ca  =  2  J  *  ,  the  total 

electric  charge  density  and  conduction  current,  are  the  stuns  of  the  partiEl 

contributions  Q  and  J*  . 

ea  •  a 

If  we  multiply  equation  (3.  18b)  by  e#/mo  and  sum  over  all  the  follow¬ 
ing  equation  emerges , 


dJ*  .  a 

TF-*  J  Tr 


o  +  — 
o  Jr 


-  I  I  —  j*  xB 


+  9. 


I  jp  .1  X  Q  (E  +  Ic  xB) 
a  ma  a  a  m«  ea  c  ° 

(3-21) 


dc  d  j 

_ o  d  _  _  •  cJ 

-ar-  +  J  -«'o- 


dt 


Spitzer  and  others  [15]  have  interpreted  the  latter  equation  as  a  genera¬ 
lized  Ohm's  law.  In  so  doing  they  have  assumed  that  the  collision  contri¬ 
bution  is  proportional  to  the  current  density  j*  when  the  state  is  not  ap¬ 
preciably  different  from  the  state  of  equilibrium,  i.-e. 


2 

a 


l  j 

c  J( 


d  t 


■  <*  sc 


a  ft 


(3.22) 


where  0  is  the  electric  conductivity.  In  the  case  that  we  neglect  time 
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variations,  quadratic  terms  in  cq,  j*  and  their  gradients,  pressure 
gradients  and  the  term  x  B  ,  equation  (3.21)  reduces  to  the  familiar 
Ohm's  law, 


j  *  =  <J(E  +  — cxB)  =  0  E*  . 

C  o 


(3.23) 


The  asterisks  denote  field  quantities  that  refer  to  a  coordinate  system 
moving  with  the  velocity  c  . 

H.  S.  Green  [8]  has  asserted  that  the  assumption  (3.22)  is  incorrect 
except  in  the  special  case  that  all  temperature  gradients  are  zero,  and  has 
proposed  another  generalization  of  Ohm's  law.  The  situation  is,  however, 
far  from  clarified,  and  we  shall  return  to  it  in  Section  19. 

Equations  (3.  18)  and  the  other  equations  in  this  section  that  are  derived 
from  them,  must  be  supplemented  with  Maxwell's  electromagnetic  equations 
(Section  3)  and  other  pertinent  equations.  But  it  is  easily  verified  that  these 
equations  do  not  constitute  a  complete,  solvable  set. 

In  the  electromagnetic  case  the  unknown  quantities  are,  according  to 
(3.28),  («‘  =  1,2,  ...  ,  »  )  , 


n 


C,  U  ,  E  ,  B  ,  n  C 
o  a  a  a 


Pa  ’  qa  ’  dT  Qa  Ca’  IT  Ue 


(3.24) 


The  number  of  scalar  quantities  thus  is  18x  +  9.  But  the  number  of  equa¬ 
tions  to  determine  the  unknowns,  i.e.  equations  (3.28),  Maxwell's  equations 
(Section  3),  the  relations  (3.26)  and  the  relation  I  Q  C  *  0  ,  as  well  as 
the  equation  of  state  a 

ill’ll  *  u«  =  |  “«*«  or  p«  *  jWu  3  • 

(3.24') 

a  total  number  of  6  x  +  15,  are,  for  all  x  =  1  ,  less  than  the  number  of 
unknowns,  and  we  need  at  least  12x  -  6  new  equations  for  a  complete, 
solvable  set  of  equations. 

Equation  (3.24')  follows  directly  from  the  definitions  of  the  respective 
quantities.  We  can  only  expect  it  to  give  the  real  pressure  for  ideal  gases. 
For  real  gases  there  is  an  additional  term  arising  from  the  interaction  of 
the  particles  [19]  .  We  have  actually  confined  ourselves  to  consider  ideal 
gases,  in  imposing  the  assumption  nD^  »  1  in  Section  8. 


It  serves  no  purpose  to  derive  new  equations  from  the  transport  equa¬ 
tion  (1.57)  by  substituting  for  $  new  molecular  properties  that  depend  on 
higher  powers  of  the  velocity.  This  method  breaks  down  because  increas¬ 
ingly  higher  moments  of  the  distribution  function  and  new  collision  contri¬ 


ve  $a 

butions  na  — j- —  are  sucessively  introduced  into  the  set  in  numbers  ex¬ 
ceeding  the  number  of  new  equations. 


It 
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19.  THE  PHENOMENOLOGICAL  RELATIONS.  OHM'S  LAW 

The  transport  phenomena  arise  when  the  gas  recedes  from  statistical 
equilibrium  as  a  result  of  disturbing  interaction  with  the  environment.  If 
this  interaction  is  of  a  transient  nature,  the  transport  phenomena  will  be 
transient.  They  will  die  away  after  a  time  equal  to  the  natural  relaxation 
time  of  the  system  tr  .  If  the  disturbances  are  constant  (time  independent), 
a  stationary  state  will  be  established  different  from  the  equilibrium  state, 
with  transport  phenomena  that  are  independent  of  time.  In  general,  condi¬ 
tions  will  vary  with  time . 

The  disturbances  may  be  represented  by  the  deviations  in  the  tempera¬ 
ture,  velocity  field,  density  and  force  field  from  the  conditions  that  are 
characteristic  of  a  gas  in  equilibrium  (equations  (2.91)  and  (2.92)).  These 
deviations  are  the  following, 


dr 


e  . 


din* 

dr 


(a  =  1,2,  ...) 


(3.25) 


where  (Section  16) 


dr 


d  In  n, 

a 

dr 


(3.26) 


Like  the  transport  vectors,  the  above  three  quantities  are  all  zero  in  the 
state  of  equilibrium  (characterized  by  Maxwellian  distribution  function),  and 
they  are  all  small  quantities  for  sufficiently  small  deviations  from  this  state. 
We  shall  hereafter  use  the  denotation  "forces"  for  the  quantities  (3.25). 

We  must  be  able  to  express  mathematically  that  the  deviation  from 
equilibrium  may  be  said  to  be  the  cause  of  transport  phenomena,  by  re¬ 
garding  all  transport  vectors  or  "fluxes"  as  functions  of  the  "forces".  A 
connection  exists  between  the  components  of  the  "fluxes"  and  "forces". 

The  reason  for  this  is  obviously  that  the  deviation  in  the  distribution  func¬ 
tion  from  the  equilibrium  distribution  is  a  function  of  the  forces . 

We  may  circumvent  the  difficulties  mentioned  in  the  preceding  section 
and  obtain  a  complete  set  of  macroscopic  equations  by  limiting  ourselves 
to  states  that  are  sufficiently  close  to  equilibrium ,  so  that  we  may  assume 
the  validity  of  a  linear  theory.  We  require  that  the  variation  in  all  macro¬ 
scopic  quantities  over  distances  of  the  order  of  magnitude  of  the  correla¬ 
tion  length  D  be  infinitesimal  [8]  .  The  distribution  functions  of  the  first 
or  higher  orders  will  then  depend  linearly  on  the  forces  (3.25). 

If  we  write  for  f#  , 


f  _  f(°)  +  f(l)  £(1)  -  £(o)  _  (1) 

~  1  a  •  1  a  ~  r, 


(3.27) 


r(o) 


where  f^w/  is  the  equilibrium  distribution  (2.73), 


,0) 


d  In  * 

d  r 


+  B. 


*  + 


we  must  have 

dr 


(3.28) 


We  require  that  f V  ^  must  not  contribute  to  the  density,  temperature 
and  mean  velocity,  so  that  the  macroscopic  parameters  na  ,  *„  and  c 

contained  in  f  '°'  are  the  real  density,  temperature  and  velocity, 
a 
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The  vectors  A  and  D  and  the  tensor  B  are  functions  of  the  velocity 
C  a  and  the  magnetic  field  B  .  They  can  be  determined  by  substitution  of 
(3.27)  in  Boltzmann's  equation  as  described  by  Chapman  and  Cowling  [4]  , 
when  their  form  is  known.  We  shall  find  the  latter  in  Section  21. 

In  a  similar  way  we  must  expect  that  the  next  approximation  to  f.  , 

12) 

i.e.  fva  ' ,  will  depend  linearly  on  terms  that  are  quadratic  in  the  forces 
and  their  gradients. 

It  follows  immediately  from  (3.28)  that  the  following  "phenomenologi¬ 
cal  relations"  must  be  valid, 


(*) 


.  2  A  •  iiii- 
dr  p  «0  dr 


(b) 

(c) 


o 


=  -2|* 


d  In  A  ' 

If- 


_iet 

dr 


(3.29) 


The  diffusion  coefficient  tensors  Aa  and  &ap  are  both  of  the  second 
order,  and  consequently  have  3^=9  components.  The  same  is  the  case 
for  the  heat  conduction  coefficient  tensors  >xa  and  )xap  •  The  viscosity 
coefficient  tensor  pa  is  of  the  fourth  order,  and  consequently  has  3^  =  81 
components . 

The  transport  vectors  QaCa  and  qa  of  course  only  depend  on  the 
i motors  ^  V1  -  and  .  ,  while  the  tensor  to  can  only  depend  on  the 

dr  dr 

tensor  e  .  The  reason  for  this  is  that  the  "cause"  and  "effect"  must  have 
the  same  symmetrical  properties  (Curie's  theorem). 

From  equation  (3.29)  and  the  relation  J*  =  e#[nftC||  it  follows  more¬ 
over  that 


J 


* 

a 


.  .  2  ^  • 

dr  J  mB  ®P  dr 


(3.30) 


The  relations  (3.29)  give  us  12  k  new  scalar  equations  for  the  macrosco¬ 
pic  quantities  (3.23)  in  addition  to  the  ones  we  already  have,  and  this  is 
sufficient  to  complete  the  set  of  equations.  (The  number  of  sufficient 
equations  was  12  h  -  6). 


By  summing  (3.29)  over  all  a  we  get  the  relations 


(•) 

a  * 

I  =  ©  , 

a  aP 

(b) 

o 

p  =  -  2  p  :  e  , 

(3.31) 

(c) 

„  _  d  In  A 

-  -  Dxl 

p  p 

dn0 
dr  ' 

where 

now 

=  I  p  ,  m 

a 

=  T  and 

a  0 

”'p  *  fm'*e 

(3.32) 
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are  the  new  coefficient  tensors . 

In  the  same  manner  we  derive  from  (3.30),  (j*  =  2  $*a  )  , 

j  *  =  -  «'  •  - ln  *  -  I  «*  •  -ini.  ,  (3.33) 

*  6r  p  0  dr 

where  the  conductivity  tensors  ® '  and  Op  are  defined  by 


Equation  (3.33)  is  the  generalized  Ohm's  law. 

In  the  special  case  that  the  gas  does  not  rotate  (-^r  CQ  =  <D )  ,  and 

we  only  have  electric  and  magnetic  forces,  equation  (3.33)  becomes  identi¬ 
cal  to  the  generalized  law  already  proposed  by  Green  [8]  ,  apart  from  the 
conductivity  in  the  latter  law  being  a  scalar. 

The  condition  is  that  we  assume  the  value  (2.  102)  for  the  chemical  po¬ 
tential,  meaning  that  the  gas  is  assumed  to  be  ideal.  This  corresponds  to 
the  assumption  made  in  Section  8  that  the  kinetic  particle  energy  is  much 

larger  than  the  interaction  energy  nD"*  »  1  .  This  assumption  is  highly 
justified  in  the  regions  of  temperature  and  density  pertinent  to  plasma 
physics.  It  corresponds  to  neglecting  the  correction  term  AFQ  (Section  9) 

which  is  identical  to  the  so-called  virial  term.  This  is  only  noticable  for 
real  gases . 

There  is,  however,  a  slight  difference  between  our  Ohm's  law  (3.33) 
and  Green's  law.  Instead  of  the  magnetic  force  Oq  x  B  contained  in 

(3.26),  Green  assumes  the  force  0&o  x  B  .  He  asserts  that  the  Hall 

flow  is  proportional  to  the  latter  quantity,  not  with  the  former.  But  when 
assuming  with  Green  that  B  =  0  in  the  equilibrium  state,  the  difference 
between  the  quantity  0Q  x  B  and  c#  Q  x  B  will  be  a  second-order  quanti¬ 
ty  because  of  the  equilibrium  condition  Oao  =  Cq  (2.91).  This  quantity 
must  be  neglected  in  a  first-order  theory. 

As  regards  Green's  criticism  of  Spitzer's  assumption  (3.22),  we  note 
that  the  equation  which  Green  identifies  with  Spitzer's  equation  (3.21)  is 
valid  for  real  gases,  while  Spitzer  implicitly  appears  to  have  assumed  that 
the  gas  is  ideal,  in  that  his  point  of  departure  is  Boltzmann's  equation 
(without  correcting  virial  term).  When  comparing  the  two  equations  it  is 
therefore  necessary  to  confine  oneself  to  ideal  gases,  especially  because 
Green's  argumentation  against  the  assumption  (3.22)  is  based  on  Chapman 
and  Cowling's  results,  which  are  only  valid  for  ideal  gases. 

For  an  ideal  gas,  where  the  collision  term  disappears  in  the  equili¬ 
brium  state,  and  thus  also  the  collision  contribution  (3.22),  it  is  evidently 
correct  to  develop  this  as  a  linear  combination  of  the  "forces"  (3.25)  when 
the  state  is  sufficiently  close  to  equilibrium.  This  is  of  course  also  the 
case  for  j*  (equation  (3.33)),  and  the  collision  term  (3.22)  may  conse¬ 
quently  in  general  be  written  as  a  linear  combination  of  j*  and  terms  in¬ 
volving  the  gradients  of  temperature  and  chemical  potential.  Only  in  the 
special  case  that  these  terms  disappear  can  Spitzer's  assumption  (3.22)  be 
correct,  and  only  provided  that  the  conductivity  is  a  tensor. 
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A  possible  criticism  must  therefore  also  be  directed  against  the  fact 
that  Spitzer  assumes  that  the  conductivity  is  a  scalar.  It  is  another  matter 
whether  there  is  really  conformance  between  the  actual  values  of  the  compo¬ 
nents  of  this  tensor  and  the  components  of  the  real  conductivity  tehsor  in 
(3.33).  This  cannot  be  established  before  the  coefficient  tensors  in  the 
phenomenological  relations  have  been  evaluated  by  direct  solution  of  Boltz¬ 
mann's  equation. 

If  equation  (3.21)  is  really  the  generalized  Ohm's  law  in  the  first  ap¬ 
proximation,  and  the  assumption  (3.22)  is  correct,  we  must,  if  we  substitute 
(3.22)  in  (3.21),  to  a  first  approximation  get  an  equation  which  is  equivalent 
with  (3.33).  The  result  course  depends  on  the  equilibrium  conditions  assu¬ 
med,  as  they  determine  which  "forces"  we  are  to  choose. 

With  the  assumption  of  his  special  equilibrium  conditions  (and  wrong 
Hall  current),  Green  believes  to  be  able  to  establish  that  for  a  binary  gas 
equation  (3.21)  cannot  be  identical  with  Ohm's  law  (3.33)  except  in  the 

d  In  X 

special  case  that  the  temperatures  gradient  — — —  is  zero. 

The  case  is  by  no  means  clarified,  however,  and  a  closer  investigation 
is  necessary. 


20.  SYMMETRY  PROPERTIES  OF  THE  MEDIUM 


When  we  are  to  determine  the  form  of  the  coefficient  tensors  A  ,  JB  , 
D  ;  W  ,  A  ,  (p  ,  in  the  preceding  section,  we  must  take  certain  assump¬ 
tions  about  the  symmetry  properties  of  the  gas. 


Position -dependent  force  fields  will  have  no  influence  on  the  symmetry 
properties,  and  will  only  be  able  to  make  the  medium  inhomogeneous.  This 
is  the  case  for  all  electric  fields,  gravitational  fields,  etc.  As  concerns 
the  only  velocity -dependent  force  entering  into  our  treatment,  the  Lorentz 


force 


1  a 

—  - (c_  xB)  ,  the  circumstances  are  different.  The  magnetic  field 

c  ma  ■ 


B  stresses  a  direction  in  the  gas,  and  all  molecules  describe  helical  trajec¬ 
tories  with  this  direction  as  an  axis .  The  magnetic  field  thus  introduces  a 
systematic  molecular  anisotropy,  which  was  not  present  in  the  absence  of 
the  field.  As  long  as  the  gyration  radius  is  sufficiently  small,  we  must 
therefore  assume  that  the  ionized  gas  has  the  same  symmetry  properties 
as  an  anisotropic  medium  with  rotational  symmetry  about  an  axis.  When 
the  magnetic  field  decreases  the  medium  will  become  gradually  isotropic. 


We  shall  investigate  the  consequences  of  this  assumption  for  the  form 
of  the  coefficient  tensor  in  the  phenomenological  relations  in  the  preceding 
section.  During  this  investigation  we  shall  need  some  results  that  will  be 
deduced  from  the  theory  of  Cartesian  tensors .  (Those  not  interested  in  the 
details  of  the  derivation  may  pass  on  at  once  to  Section  22,  just  accepting 
the  validity  of  the  results  of  Section  21). 


M 
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2  1 .  CARTESIAN  TENSORS 


A  Cartesian  tensor  T  of  the  p-th  order  is  defined  in  a  3- 

r lr2  ‘  ‘  rp 

dimensional  Euclidean  space  as  a  set  of  3?  quantities  (numbers,  coordinate 
functions,  etc.)  that  by  an  orthogonal  transformation  of  coordinates, 


x*i  =  aikxk  +  bi  ' 


aikajk  =  5ik  • 
is  transformed  in  accordance  with  the  rule 


aik  I  =  +  1  ’ 


(3.35) 


T  ili2  -  •  lp  ailriai2r2  &iprp  lr2  •  ‘  •  rp 


(3.36) 


Summation  is  everywhere  performed  from  1  to  3  over  indices  that  occur 
twice  in  the  same  term.  The  old  and  the  new  coordinates  are  denoted  by 
the  vectors  x^  and  x1^  .  A  vector  is  a  tensor  of  the  1-st  order,  and  a 
scalar  is  a  tensor  of  the  0-th  order.  The  transform  of  a  quantity  is  deno¬ 
ted  by  an  asterisk. 


An  orthogonal  transformation  (3.35)  is  always  composed  of  a  pure 
translation  (bj  =  Const.)  and  a  pure  rotation  (a^3^’  a^  =  Const.)  of  the 

coordinate  system.  Rotation  through  an  angle  ♦  about  an  axis  with  direc¬ 
tion  vector  1^  can,  according  to  Spain  [20]  ,  be  represented  by  the  trans¬ 
formation  matrix, 


lk 


=  cos  6-k  +  (1  -  cos  O^lfc  +  8in'i’eijk1j 


(3.37) 


Here 


123 


=  e 


312 


=  c 


231 


+  1, 


:32  1 


all  the  other 


=  0.  It  can  be  shown  that 


*213  ” 
&ik  and 


132 

'ijk 


=  -  1 ,  and 
are  tensors 


og  the  second  and  third  order.  If  the  angle  of  rotation  4  is  infi  litesimal, 
+  =  6  4  >  the  relation  (3.37)  turns  into 


aik 


*ik  + 


**Cijklj 


=  6ik  +  Sik 


(3.38) 


As  the  inner  product  of  the  two  tensors  is  a  new  tensor,  then  a.^  and  S^ 
must  also  be  tensors,  according  to  (3.37)  and  (3.38). 


21.1  TRANSFORMATION  OF  LINEAR  TENSOR  RELATIONS 

Suppose  that  a  linear,  homogeneous  relation  exists  between  the  Carte¬ 
sian  tensor  T_  _  and  another  tensor  U_  _  (of  the  order  q), 

.  (3.39) 


r . .  . .  r 
1  P 


■!••• sq 


=  K 


u 

rp»  *  1  •  •  •  Sq 


where  T,  a  and  U  are  functions  of  the  coordinates  x^  and  the  time  t. 

According  to  the  quotient  theorem  for  tensors ,  the  3^*^  coefficients 
are  components  of  a  tensor  of  the  order  p+q  .  (The  quotient  theorem 
states  that  a  set  of  3r  quantities  are  components  of  a  tensor  of  the  order 
r  if  the  inner  product  of  these  quantities  with  a  tensor  with  arbitrary  com¬ 
ponents  is  also  a  tensor.)  The  transformation  rule  for  the  coefficient  ten¬ 
sor  is  therefore 
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=  ai .  r  ‘  •  •  ai  r 
11  p  p 


a.  ...  a.  k 

Jlal  Vq  r,-‘r 


r 


P.8r--8q  • 


(3.40) 


The  dependence  (3.39)  must  therefore  also  exist  between  the  transformed 
quantities.  The  relation  (3.39)  is  invariant  with  respect  to  coordinate 
transformations . 

For  certain  transformations  a^,  dependent  on  the  symmetry  pro¬ 
perties  of  the  medium,  the  tensor  K  will  be  transformed  into 

rl . . . sq 

itself.  If  the  medium,  for  instance,  is  isotropic,  this  will  be  the  case  for 
all  rotations  of  the  coordinate  system;  if  the  medium  is  rotationally  sym¬ 
metrical,  for  all  rotations  about  the  axis  of  symmetry,  etc.  For  a  homo¬ 
geneous  medium  all  coefficients  x  will  be  constant  (invariance 

r  1  ‘  '  *q 

with  respect  to  a  pure  translation  of  the  coordinate  system). 

If  we  know  the  symmetry  properties  of  the  medium,  we  thus  know  the 

transformation  tensors  a.,  for  which  the  relation 

lk 


*i  •  •  *  Vh 


(3.41) 


is  satisfied.  Together,  expressions  (3.40)  and  (3.41)  then  give  a  certain 

number,  namely  3^+q  •  S,  of  equations  to  determine  the  quantities 
x  The  quantity  S  is  the  number  of  transformations  for  which 

r 1  ‘  ‘  8q 

(3.41)  is  satisfied.  Only  a  certain  number  of  these  equations  will  be  in¬ 
dependent  of  each  other,  so  that  we  cannot  expect  to  be  able  to  determine 

all  the  3**+t*  coefficients  completely.  What  we  can  determine  by  this  pro¬ 
cedure  is  the  form  of  the  coefficient  tensor  which  corresponds  to  the  sym¬ 
metry  properties  of  the  medium. 

We  shall  now  carry  out  this  line  of  reasoning  in  detail  for  the  tensor 
relations 


Ti  * 


x.  U  , 

ir  r 1 


Tu  1 


x..  U 
ijrs  rs 


(3.42) 


under  the  assumption  of  a  rotationally  symmetrical  medium.  The  quanti¬ 
ties  T^j  and  Urg  are  both  assumed  to  be  symmetrical. 

For  simplicity  we  place  the  x^ -axis  of  our  coordinate  system  along 
the  axis  of  symmetry.  The  requirement  of  rotational  symmetry  then 
means  that  the  relation  (3.41)  must  be  satisfied  for  arbitrary  rotations 
of  the  coordinate  system  about  the  X3-axis.  As  this  axis  has  the  direction 
vector  1^  =  the  transformation  tensor  must,  according  to  (3.37), 

in  this  case  have  the  form 


*ik  =  cos  ♦  $ik  +  (1  -  cos  +  )  6i3  4k3  +  sin  ♦  *i3k 


(3.43) 


or,  on  matrix  form, 
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Caik3 


cos  «(*  -Bin  $  0 

sin  t  cos  <|*  0 

0  0  1 


(3.44) 


The  requirement  of  invariance  (3.41)  for  *x, 

lT 


is  satisfied  for  arbi 


trary  rotations  ♦  about  the  x^-axis  if  it  is  satisfied  for  the  three  rota¬ 
tions  ♦  =  *,  ♦  =  x/Z  and  4»  =  64>  with  the  transformation  matrices , 
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(3.45) 


21.2  THE  RELATION  T.  =  x.  U 

1  llr  T 

The  coefficient  tensor  x.,  has  the  matrix 

lk 


Tl 

H12 

*13 

21 

*22 

*23 

31 

*32 

*33 

(3.46) 


In  this  case  the  invariance  requirement  (3.41)  gives 

i..  or  (a.  a.  -  6.  &.  )n 
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In  this  equation  we  successively  substitute  the  three  sets  (3.45)  for  (a.,  ) 
and  find  the  corresponding  conditions  for  the  components  xrg  . 


The  requirement  of  invariance  with  respect  to  the  rotation  f  «  g 
gives  the  conditions  xj3  -  x^ 

duces  to 
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The  requirements  of  invariance  with  respect  to  the  rotation  +  =  x/z 
further  gives  the  new  conditions  =  ~  *\2’  *22  =  *il»  and  [Hik] 

is  further  reduced  to 
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The  requirement  of  invariance  with  respect  to  the  rotation  f  =  64  means 
that  equation  (3 . 47)  is  satisfied  with  a^  =  6^k  +  S^, 

We  have,  in  other  words,  the  relation 
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when  quadratic  terms  in  S.,  are  neglected.  These  equations,  however, 
give  no  new  conditions,  ana  (3.49)  thus  is  the  most  general  rotationally 
symmetric  tensor  of  the  second  order. 


If  we  so  wish,  we  can  write  the  relation  T^ 
vector -dyadic  form  as  T  =  (|* .  U  ,  ||*  =  *>.  x^ 

unit  vectors  in  the  three  directions  of  the  axes, 
pression  for  T  can  be  transformed  to 


=  x  .  U  on  the  familiar 
ek  where  e.  denotes  the 

If  we  use  (3.49),  the  ex- 


T  =  *n  +  »12(0  x  e3)+  x33Uf|  =  ix  •  D  ,  (3.51) 


where  0^=  Ui®i  +  U2«2  *nd  ®||  =  U3  ®3  are  the  components  of  the 
vector  U  normal  to  and  parallel  with  the  axis  of  symmetry  (  •  3) . 

Finally  we  shall  only  note  that  the  most  general  iaotropto  tensor  of 
the  second  order  has  the  form 
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(3.52) 


This  can  be  easily  verified  by  requiring  invariance  even  for  arbitrary  rota¬ 
tions  about  the  x  j  -  and  X2~axes. 


21.3  THE  RELATION  T..  =  x ..  U 

ij  ijrs  rs 

4 

In  this  case  the  coefficient  tensor  has  3  =81  components,  which  may 

be  arranged  on  matrix  form  according  to  the  following  convention, 


IT 
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If  we  now  write  Tjj  and  U^j  as  column  matrices  [T^]  and  [U^]  where 

the  components  are  arranged  above  each  other  according  to  the  following 
convention , 
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(3.54) 


and  correspondingly  for  U..  ,  the  relation  T..  =  »,.  U  can  be  written 
on  matrix  form  ..  1J  *J  l,r*  ” 


tTij)  “  '  IDr.l  ' 

Because  of  the  symmetry  properties 
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(3.56) 


only  36  of  the  81  coefficients  *. ,  are  really  independent.  It  follows  from 
(3.56)  that  1JW 


*..  U  =  *..  U  =  U  =  U 
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(3.57) 


and  as  the  components  may  be  arbitrarily  chosen,  we  must  have  that 


*ijrs  "  *ijsr  =  *jirs  "  *jisr  ' 


(3.58) 


Thus  H^jra  i*  symmetrical  with  respect  to  the  indices  i  and  j  ,  and  the 
indices  r  and  s. 

The  independent  coefficients  are  the  6x6  =  36  in  the  upper  left-hand 
corner  of  (3 . 53).  The  coefficients  *3231’  *2313’  *3213  are*  *or  instance, 

all  equal  to  N2331'  note  that  if  Up#  is  anti-symmetrical,  then  also 
x^jrt  is  anti -symmetrical  in  the  indices  r  and  s  , 
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X..  =  X.. 

ijrs  jirs 


*ijsr  =  "  *jisr 


(3.59) 


Instead  of  writing  the  whole  matrix  (3. 53)  we  will  henceforth  only  write 
down  the  36  independent  elements  in  the  upper  left-hand  corner. 

The  invariance  requirement  now  gives  us 
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(3.60) 


We  again  substitute  successively  the  three  transformations  matrices 
(3.45)  representing  the  rotations  about  the  x^-axis,  and  obtain  correspond¬ 
ing  conditions  for  the  coefficients  x^j^  ■ 

The  requirement  of  invariance  for  the  rotation  #  =  t  gives  the  condi¬ 
tion 


when  index  3  occurs  an 
uneven  number  of  times. 


(3.61) 


This  means  that  the  6  X  6  -  matrix  for  the  independent  quantities  x,...  get 
the  following  reduced  form,  ^ 
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The  number  of  independent  coefficients  is  thus  reduced  from  36  to  20. 

The  requirement  of  invariance  with  respect  to  the  rotation  4  =  s/2 
further  gives 
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The  number  of  20  independent  coefficients  further  reduces  to  10,  and 
(3.62)  gets  the  form 
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The  invariance  requirement  for  the  rotation  4  =  &♦  with  the  trans- 

(3.65) 


formation  tensor  aik  =  6^  +  s^  >  s^  =  6**^3^  gives 
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which  leads  to  the  relations 
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The  final  form  of  the  rotational  symmetric  coefficient  matrix  thus  becomes 
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1112  1112  w 

where  only  8  of  the  coefficients  are  independent. 

For  the  sake  of  comparison  we  shall  also  give  the  result  for  an  tao- 
tropio  medium  (invariance  with  respect  to  rotation  about  all  three  coordi¬ 
nate  axes) : 


ijkl 

or  on  matrix  form, 
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This  is  the  most  general  isotropic  tensor  of  the  fourth  order,  which  thus 
has  three  independent  coefficients. 


22 .  PHENOMENOLOGICAL  RELATIONS  FOR  ROTATIONALLY 
SYMMETRIC  MEDIUM 

In  the  following  we  shall  assume  without  further  verification  that  the 
vectors  and  tensors  we  have  worked  with  in  the  preceding  sections  are  real¬ 
ly  Cartesian  tensors,  i.e.  they  are  transformed  in  accordance  with  (3.36). 
We  further  assume  the  symmetry  properties  of  the  medium  adopted  in  Sec¬ 
tion  20,  and  that  all  tensors  we  have  to  do  with  are  rotationally  symmetri¬ 
cal  (or,  in  the  non -magnetic  case,  isotropic). 

Let  us  first  consider  the  correction  (3.28)  to  the  distribution  function. 
We  suppose  that  the  vectors  A  end  D  are  linear  functions  of  the  thermal 
velocity  C  ,  and  that  the  tensor  E  is  a  linear  function  of  the  symmetric 

and  non-divergent  tensor  C°C  =  C  C  ■  |  C^  I ,  i.e.  that 

A=«(C)C,  £  =  P(C)  :  C°C  ,  D  =  v  (C)  ■  C  ,  (3.70) 

(«  ,  and  jf  are  isotropic  tensors  for  an  isotropic  medium  (  B  =  0  )).  As 
the  two  second-order  tensors  tt  and  g  and  the  fourth-order  tensor  (0  are 
rotationally  symmetrical,  and  the  former  two  consequently  have  the  form 
(3.49),  while  fi  has  the  form  (3.67),  we  have,  according  to  equation  (3.51), 
that  A  and  0  are  linear  combinations  of  the  vectors 

C±  C  x  m°  C|(  (3,71) 

(«s°  is  the  unit  vector  in  the  direction  of  the  magnetic  field),  or,  which  is 
equivalent ,  of  the  vectors 

C  C  x  •«  (C  x  u°)  X  «s°  .  (3.72) 

A  closer  investigation,  the  details  of  which  we  shall  omit  here,  shows  that 
when  0  has  the  form  (3.67),  then  B  can  be  written  as  a  linear  combina¬ 
tion  of  the  six  independent,  symmetric  and  non-divergent  tensors  formed  by 
combination  of  the  vectors  (3.71)  or  (3.72). 

If  the  vectors  (3. 72)  are  denoted  by  lj,  ^  and  1^,  respectively,  the 
six  tensors  involved  are  given  by  the  expressions 

i(D#  ttfi  +  Ip  ■,),  «.*  *  1.2.3  .  (3.73) 

All  coefficients  in  the  relations  (3.70)  are  in  general  functions  of  the 
scalar  value  C  of  C  . 

In  a  similar  way  we  assume  the  tensors  »  ,  A ,  0  and  pi  in  the 
phenomenological  relations  (3.31)  and  (3.33)  to  be  rotationally  symmetrical, 
i.e.  that  the  three  first  have  matrices  of  the  form 


U 
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while  the  last  tensor  can  be  represented  by  the  matrix 
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(3.75) 


Digressing  briefly,  we  note  that  (3.74)  also  is  the  form  of  the  dielectric 
tensor  c  in  the  phenomenological  relation  D  s  £.  .  E  which  is  generally 
assumed  to  be  valid  in  the  study  of  wave  propagation  in  ionized  gases.  (See 
for  instance  [3]  .)  This  appears  to  support  our  assumption  about  the  sym¬ 
metry  properties  of  the  medium. 

Between  the  coefficients  in  the  phenomenological  equations  there  must 
exist  relations  that  can  be  derived  by  the  use  of  the  so-called  principle  of 
microscopic  reversibility  [Zl]  .  These  relations  are  the  Onsager  rela¬ 
tions  . 

We  shall  not  enter  further  into  this ,  especially  because  the  conditions 
do  not  appear  to  be  quite  clarified  yet,  but  only  note  that  we  must,  among 
other  things,  have 


*ij(B)  =  *ji(-B)  ,  *..(B)  =  X^-B) 

V*)  =  V-B>-  >m{B)  s  W-B>  • 


(3.76) 


This  means  that  the  diagonal  elements  and  all  symmetric  elements  in 
(3.74)  and  (3.75)  must  be  even  functions  of  the  magnetic  field  B  ,  while 
the  anti -symmetric  elements  must  be  uneven  functions  of  the  field.  These 
conclusions  are  in  accord  with  the  results  that  may  be  derived  by  direct 
solution  of  Boltzmann's  equation  (Section  25). 

If  w  ,  B ,  V  and  yt  had  been  isotropic  tensors ,  the  different  flux 
vectors,  diffusions,  electric  current,  heat  of  flow  and  flow  of  momentum 
would  have  been  directly  proportional  to  the  corresponding  forces.  This 
is  usually  the  case.  In  our  case,  on  the  other  hand,  this  is  no  longer  so. 
Fluxes  are  no  longer  parallel  to  the  forces.  A  temperature  gradient 
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—  and  a  force  will,  for  instance,  give  origin  to  diffusion,  elec¬ 

tric  current  and  heat  flow  in  the  three  directions 
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or  in  the  directions 
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At  the  same  time  the  usual  viscosity  hypothesis  [22]  which  assumes  that 
the  stress  tensor  is  directly  proportional  to  the  deformation  tensor 

(£  =  -2|i  c)  breaks  down.  This  assumption  is  only  correct  for  an  iso¬ 
tropic  medium  (  B  =  0  /  and  must  be  supplanted  by  the  more  general  hypo- 
p  =  >2  (i  :  e  where  p  now  has  the  form  (2. 75). 

For  the  time  being  the  phenomenological  coefficient  tensors  (  M ,  B, 
«  and  (p  )  are  unknown  parameters  in  a  set  of  macroscopic  equations.  In 
Sections  24  and  25  we  shall  compute  these  quantities  directly  with  the  aid 
of  Boltzmann's  equations. 


23.  IRREVERSIBLE  PROCESSES 


The  macroscopic  equation  for  the  internal  energy,  defined  as  average 
kinetic  energy  per  unit  volume,  is  (3.  19c), 


°o  +  dr 
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q  +  P  : 


-Z  4  C  f  *  o  .  (3.77) 
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(F#  now  contains  the  magnetic  force  CQX  )  .  This  equation  is  in  accord 

with  the  energy  equation  given  in  [8]  ,  but  is  derived  in  another  and  more 
simple  manner. 

The  usual  macroscopic  thermodynamic  entropy  equation  for  quasi¬ 
static  processes, 
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(3.78) 
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does  not  yet,  like  (3.77),  have  a  microscopic  justification.  In  this  (Gibb's) 


equation,  S,  U  and  denote  entropy,  internal  energy  and  chemical 

potential  per  unit  mass,  respectively.  For  an  ideal  gas  is  now  given 

by  (2.  101).  We  introduce  instead  the  corresponding  quantities  per  unit  vo¬ 
lume,  and  with  the  aid  of  the  equation  of  continuity  (3.  18a),  equation  (3.78) 
gets  the  form 
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(3.79) 

Combination  of  (3.78)  with  (3.77)  gives  us  the  following  expression  for  the 
entropy  change, 
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If  we  introduce  the  entropy  flow 


q„  =  <1  "  l*an*Ca 

ft 


(3.80) 

(3.81) 


and  the  generalized  chemical  potential  (2.  107) 
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(3.82) 


where  4a  t^ie  total  potential  for  all  force  fields  (gravitational  fields, 
electromagnetic  fields,  centrifugal  fields,  etc.),  equation  (3.80)  can  be 
written 
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(3.83) 


The  contribution  to  this  equation  from  all  forces  that  are  proportion^}  to  the 
mass  (  4a  independent  of  a  )  is  zero  because  of  the  relation  Z  Q#Ca  s 

a 

This  holds,  for  instance,  for  gravitational  forces  and  centrifugal  forces. 

If  we  reintroduce  entropy  per  unit  mass,  (3.83)  gets  the  form 
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In  this  equation  of  continuity  the  right-hand  side  must  evidently  be  interpre¬ 
ted  as  the  source  density  of  the  entropy.  If  we  introduce  for  the  "fluxes"  q, 
o  _ 

p  and  QgCa  the  expressions  in  Section  1 9 »  this  source  density  will  become 

a  quadratic  form  in  the  "forces"  —  — —  •  ,  «  and  .  We  must  require 

dr  dr 

that  this  form  is  a  positive  definit.  It  can  then  be  interpreted  as  the  irrever¬ 
sibly  dissipated  amount  of  heat  per  unit  volume  and  time.  The  expressions 
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give  the  respective  thermal  dissipation,  viscous  dissipation  and  dissipation 
because  of  gradients  in  the  generalised  chemical  potential.  The  last  term 
will  contain  the  Joule  dissipation  •  E  *  (3.23).  As  we  have  seen  in 

Section  22  the  connection  between  "fluxes"  and  "forces"  is  strongly  modified 
by  the  presence  of  a  magnetic  field.  The  same  will  consequently  be  the  case 
for  the  dissipation. 

As  we  have  assumed  the  validity  of  (3.78)  during  the  derivation  of  (3.84), 
the  latter  equation  has  the  character  of  a  macroscopic  postulate.  The  vali¬ 
dity  of  this  postulate  must  be  more  closely  investigated  by  comparison  of 
(3.84)  with  the  equation  of  continuity  (2.  14), 
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which  is  derived  on  a  microscopic  base.  We  multiply  (3.86)  by  -kT,  sum 
over  all  a  and  get  (S  =  -kH) , 
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For  quasistatic  processes  it  can  be  easily 


shown  that  the  entropy  flow 
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is  identical  to  (3.81).  We  have,  namely, 

6qs  =  -  2  kT  /Ca(l  +  lnfo)4fadce 


(3.89) 


If  we  write  6 f  =  f  -  ity  =  fi°^  f-  and  8®  =  q_  -  q^°^  and  take 

a  a  Sz.\  *  ■  ■ 

into  consideration  that  '  =  0,  (3. 89)  becomes  approximative^ , 


<i.  -  -  *  m  jc  0  + 1»  do, 


-  IkTj  0.0  *  to  do. 


If 
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I  nflkTCfl  (l  +  lnf^)  =  q-t(.a°naCa 


a  “  “  a 

We  can  therefore  write  equation  (3.87)  as 

T(4£  +  si.  c)+t1.  (I  q.)  =  T  a 

dt  dr  °  dr  t 

or,  if  we  introduce  the  entropy  per  unit  mass, 


(3.90) 


(3.91) 


QT§  +  t£- «s>  =  «Td>  0 


(3.92) 


The  entropy  source  density 


QO 


*cS 

dt 


-k  I  /  In  f  (/cf“ 

«P 


dt 


)  do 

P  « 


(3.93) 


is,  as  shown  in  Section  12  always  positive.  It  can  be  evaluated  microscopi¬ 
cally  by  substitution  of  Landau's  collision  term  in  equation  (3.93)  and  by 

writing  fa  =  f^°^  (1  +f#).  The  form  for  fa  is  (3.28).  The  coefficients 

in  this  relation  must  be  determined  with  the  aid  of  Boltzmann's  equation. 

It  is  easily  seen  that  qb  becomes  a  quadratic  form  in  the  "forces". 

Equation  (3.92)  is  completely  identical  to  (3.84)  if  this  quadratic  form 
is  the  same  as  the  right-hand  side  of  (3.84), 


Q  T  a  =  -  q  .  d-ln  -•  -  p  :  e  -  I  n  C 

^s  dr  *  ® 


(3.94) 


o  _ 

The  quantities  b  ,  q  ^  ,  jp  and  naCa  can  all  be  evaluated  when  fa  is 

determined.  By  substitution  in  equation  (3.94)  we  can  easily  verify  whether 
this  relation  is  correct  or  not.  As  far  as  we  know  the  validity  of  this  rela¬ 
tion  has  not  been  established  on  a  microscopic  basis  by  any  authors. 
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PART  IV.  SOLUTION  OF  BOLTZMANN'S  EQUATION 


24.  THE  METHOD  OF  SUCCESSIVE  APPROXIMATIONS 


According  to  whether  the  relaxation  time  for  approximation  to  statis¬ 
tical  equilibrium  for  our  system  is  essentially  less  than,  equal  to  or  greater 
than  the  characteristic  times  t  for  other  physical  effects,  we  must  distin¬ 
guish  between  three  different  cases: 

(a)  «  t  .  Collisions  are  of  dominant  significance,  the  state  is  close  to 

statistical  equilibrium. 


(b)  t  *•  t  Collisions  are  of  the  same  importance  as  other  physical  ef¬ 
fects  . 


(c)  t  »  t  .  Collisions  are  of  subordinate  importance.  The  state  is  far 
from  statistical  equilibrium. 

For  an  electron -proton  plasma  we  can  put  Tr  equal  to  the  equalization 
time  Tj  for  momentum  and  energy  among  the  electrons  (2.48).  By  substi¬ 
tution  of  numerical  values,  the  conditions  (a),  (b)  and  (c)  become  the  follow¬ 
ing. 


(4.1b) 

c 


The  case  (a)  is  only  had  when  the  density  is  sufficiently  large  and  the  tempe¬ 
rature  sufficiently  small.  In  case  (c)  on  the  other  hand,  the  density  will  be 
very  small,  and  the  temperature  large. 

In  case  (b)  the  state  of  the  system  will  continue  near  statistical  equili¬ 
brium  if  t  is  the  gyration  time  for  the  particles  in  a  magnetic  field.  A 
magnetic  field  will,  namely,  not  have  any  effect  on  a  Maxwellian  velocity 
distribution.  The  condition  is  that  the  thermal  velocity  Cg  is  much  larger 
than  the  mean  velocity  cQ  so  that  we  can  write  approximative^  cft  s  Ca 

For  isotropic  velocity  distribution  ffl  =  fg  (C#  )  the  magnetic  force  term 
4f 

(c  x*»  )  •  — —  in  Boltzmann's  equation  (1.56)  will  vanish,  and  the  mag- 
°  a 

netic  field  will  only  be  of  importance  for  the  higher  approximation  to  the 
distribution  function. 

We  have  previously  seen  that  we  can  get  a  complete,  solvable  macro¬ 
scopic  set  of  equations  when  we  confine  ourselves  to  the  treatment  of  the 
cases  (a)  and  (b).  To  achieve  this  we  must  introduce  the  so-called  pheno¬ 
menological  relations  that  contain  the  unknown  parameters,  the  coefficient 
tensors.  In  case  (c)  we  can  generally  not  derive  any  complete  macroscopic 
set  of  equations,  and  we  are  restricted  to  a  microscopic  treatment. 

The  alternative  microscopic  treatment  consists  in  a  direct  solution  of 
Boltzmann's  equation.  As  this  equation  is  extremely  complicated  in  the 
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general  case  and  is  impossible  to  solve  in  an  exact  manner,  we  must  resort 
to  approximative  methods  of  solution.  The  most  usual  method  is  that  of  suc¬ 
cessive  approximations. 

This  method  can  be  described  schematically  as  a  convenient  subdivision 
of  Boltzmann's  equation  into  two  or  more  (possibly  an  infinite  number  of) 
terms,  of  which  the  mutual  ratios  become  the  expansion  parameter  in  a 
series  development  of  the  distribution  function  fQ  : 


By  substitution  of  this  expression  in  Boltzmann's  equation,  the  1st.  2nd,  etc. 
approximations  to  fB  are  determined.  Which  quantity  to  choose  for  the  ex¬ 
pansion  parameter  naturally  depends  on  the  nature  of  actual  physical  problem 
and  the  assumed  conditions.  This,  that  different  states  for  the  system  require 
different  procedures,  is  of  course  a  great  weakness  of  this  method. 

In  case  (a)  and  (b),  when  inter  molecular  interactions  (and  in  case  (b)  the 
effect  of  the  magnetic  field)  are  the  dominant  physical  effects ,  the  method  of 
solution  of  Chapman  and  Enskog  may  be  applied.  This  method  is  described 
in  detail  by  Chapman  and  Cowling  [4]  .  In  case  (c)  the  procedure  will  depend 
upon  the  different  physical  conditions .  Chew,  Low  and  Goldberger  [23]  as¬ 
sume,  for  instance,  that  the  dominant  physical  effect  is  due  to  electric  and 
magnetic  fields,  and  completely  disregard  collisions.  This  assumption  is 
taken  as  the  point  of  departure  for  a  determination  of  the  distribution  function 
by  successive  approximations. 

In  the  following  we  shall  confine  ourselves  to  cases  (a)  and  (b)  to  make 
comparisons  with  the  macroscopic  theory  developed  in  earlier  sections  for 
near -equilibrium  states.  We  shall  not  describe  the  method  of  solution  of 
Chapman  and  Enskog,  but  only  cite  the  results  it  leads  to  in  the  1st  and  2nd 
approximation.  The  manner  of  presentation  is  essentially  that  of  Chapman 
and  Cowling.  Along  with  these  authors  we  shall  restrict  ourselves  to  the 
treatment  of  gases  with  one  or  two  components ,  and  for  the  sake  of  coherence 
include  both  the  non -magnetic  and  the  magnetic  cases. 

The  method  of  Chapman  and  Enskog  is  based  on  Boltzmann's  equation 
with  the  short-range  collision  integral  (1.26).  But  we  have  shown  earlier 
that  the  correct  form  of  the  collision  integral  for  ionized  gases  is  Landau's 
form  (1.27).  Landau's  collision  integral  does,  however,  like  (1.26)  possess 
exactly  the  formal  properties  essential  to  the  application  of  Chapman  and 
Enskog 's  method.  It  can  be  split  up  in  the  same  manner  as  (1.26);  it  is, 
like  (1.26),  zero  for  Maxwellian  velocity  distribution,  and  has  the  same 
summation  invariants  (1,  mo  ,  Jmc2  ).  We  must  therefore  assume  that 
we  can  use  Landau's  integral,  and  still  use  Chapman  and  Enskog 's  method 
of  solution. 

Considering,  however,  that  a  critical  investigation  of  the  above  is  mis¬ 
sing,  we  shall  use  the  familiar  integral  (1.26).  If  our  assumption  that  the 
two  collision  terms  will  lead  to  formally  identical  results  is  correct,  then 
the  choice  of  collision  integral  is  unimportant  as  long  as  we  are  only  inter¬ 
ested  in  general  formulas,  and  not  in  numerical  results. 

For  the  sake  of  clarity  we  write  down  Boltzmann 's  equation  with  the 
short-range  collision  integral  (1.26), 


dt 


4f. 


+  C, 


dr 


(F  +  o  x  « 
'a  o  a 


d-f£) 

dt  ' 


of 


dC_ 


+  (Ca  x  wa  ) 


df 

_ a 
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df. 


4  f 
c  a 


4Ca  ’’  dr  °o  =  '  p  Sa^f#,f0*  "  dt 
where  Sftp  i»  given  by 

S«pMp>  =  -JJ  (**  -£efp)k«p^  dcp 


(4.2) 


(4.3) 


In  the  case  (a)  we  thus  assume  that  the  collision  term  is  dominant  in 
(4.2)  (and  write  for  simplicity  B  =  0  ),  while  in  case  (b)  we  assume  that 

the  collision  term  and  the  magnetic  force  field  (Cftxwa)  •  -  -  are  the 
dominant  terms.  a 


24.1  CASE  (a).  SIMPLE  GAS 


To  determine  the  1st  approximation  f^°^  to  the  distribution  function 

*  (<r> 
o 


we  get 

S((<°>,  f<°>)  =  0  , 

which  leads  to  Maxwellian  velocity  distribution 


f(°) 


v.  ,  mc2- 


2kT 


(4.4) 


(4.5) 


The  parameters  n,  cQ  and  T  can  be  interpreted  as  the  respective  real 
density,  mean  velocity  and  temperature  if  we  require  that  the  relations 

/  4  f(r)do  *  0,  r  •  1.2  ...  (4.5a) 

* 

be  satisfied  for  each  of  the  three  summation  variants  4  *  (1,  mo  ,  jmc  ). 
To  determine  the  2nd  approximation  *  f(°)  f  ^  we  get  the  equa¬ 


tion 


V 


(o) 


+  c 


4f<°>  .  doco  „  4f<°>  4f<°>  ..  4 


+  (F 


-) 


dt  4r  dt  '  4C 

+  S(f^°^,  f(°)f(1))  +  s  (^f^\  f^)  »  0 


„  C  :  —  c  + 
dC  dr  ° 


(4-6) 


d  4  4 

Here  =  rr  +  c  •  —  .  The  operator  —  has  a  fixed  form  for  each 
dt  4t  o  4r  r  dt 

of  the  quantities  (n,  o  ,  T)  it  works  on  [4]  .  Equation  (4.6)  is  linear 


In  accord- 


both  in  the  "forces"  tfe-X  and  •  ,  and  moreover  in  e  ^ 

dr 

ance  with  Section  19.  f(*)  must  therefore  be  a  linear  combination  of  the 
forces 
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(4.7) 


9 


(1)  = 


A- 


din  T  + 
dr 


E  :  «  . 


To  determine  A  and  £  we  insert  (4.7)  into  (4.6)  and  the  coefficients  be- 
d  In  T 

fore  the  forces  — —  and  «  are  put  equal  to  zero.  In  this  way  we  get 
o  r 

separate  equations  for  A  and  £  .  According  to  these  equations,  A  and  IB 
have  the  forms  that  might  be  expected  for  an  isotropic  medium  (Section  22), 


A  =  «  (C)  C  ,  B  =  p(C)C°C 


(4.8) 


Next,  the  coefficients  a  and  0  are  expressed  as  converging  series  of 
Sonine  polynomials  (generalized  Laguerre  polynomials), 

«(C)  =  IarSr/a(tfZ).  8  (C)  =  £  brs/*  (tf2)  .  (4.9) 

1/ 

Here  the  dimension-free  quantity  %  -  (— — )  C  is  introduced  to  re¬ 

place  C 

We  finally  substitute  these  series  in  the  equations  for  A  and  £, 
multiply  the  equations  by  (&)  and  respectively,  and 

integrate  over  velocity  space.  By  using  the  orthogonality  conditions  for 
Sonine  polynomials ,  we  thereby  get  an  infinite ,  linear  set  of  equations  for 
the  coefficients  a  and  another  set  for  b^.  These  sets  can  be  solved  ap¬ 
proximately  with  the  desired  accuracy  by  breaking  off  the  series  after  a 
given  number  of  terms. 


24.2  CASE  (a).  BINARY  GAS 


To  determine  the  1st  approximation  f^°^  to  f#  = 
we  get  in  a  similar  way  as  for  a  simple  gas , 


f(r) 


I  f'*'  («=  1.2) 

o  • 


*.  S«P<W  a  0 


8=1 


a  =  1.2 


(4.  10) 


The  solution  of  these  two  equations  is  again  Maxwellian  distribution  with  the 
same  temperature  and  average  velocity  for  the  two  gases, 


*\2«kT  ) 


V. 


exp 


•(- 


m  C 

a  a 

2kT 


(4.11) 


Here  n#  oQ  and  T  are  as  before  identical  with  the  actual  physical  quanti¬ 
ties,  on  the  assumption  that  the  higher  approximations  f(r),  r  ?  1  do  not 

a  ^  (1) 

give  any  contribution.  The  equations  for  the  next  approximation  fv#  7  = 

f^)  f  are  tzinlogous  with  the  corresponding  equation  (4.6)  for  a  simple 
gas,  and  the  solution  has  the  form 
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•  0)  =  A  •  d  1-T  +  B  :  e  +  D  •  d 
Ta  a  dr  a  o  12 


(3.12) 


The  vector  dl9  is  defined  by 


‘12 

m 


12 


1 

d  ~n~  _  nln2^m2~ml^  d  In  p  _  Q1Q2  _  p 
dr  «Q  dr  pQ  1  2  ‘ 


(4.13) 


By  some  calculation  we  can  show  that  (4.  12)  is  completely  equivalent  with 
the  expression  (3.28).  We  can,  namely,  write 


■12 


where 


Q1Q2 

PQ 


1 

m. 


(di  -  d2)  , 


***  -  (u°  5  kTl  d  lnT 

—  -  -  ?kT)  -j— 


(4.  14) 


(4.15) 


As  before  and  pa  are  the  respective  ordinary  and  generalized  chemi¬ 
cal  potentials  (2.  101)  and  (2.  107).  The  last  term  in  equation  (4.  12)  can 
therefore  be  transformed  to 


2 

I  D 
P  =  1 


dp 


D«1  "  -D.2 


Ql92 

i  ■  ■  D 

PQ  « 


(4.  16) 
d  In  T 

dr 


If  we  include  the  part  of  this  term  which  is  explicitly  dependent  on 

in  the  first  term  of  (4.  12),  the  expression  becomes  precisely  equal  to 
equation  (3.28). 

In  the  same  way  as  for  a  simple  gas  we  further  find  that  the  form  of  A  , 
S  and  D  is  the  one  to  be  expected  for  an  isotropic  medium 

A  =  «(C)C,  D  =  y(C)C,  B  =  8(C)C°C  .  (4.17) 

The  coefficients  a  ,  0  and  y  are  evaluated  in  the  same  way  as  previously. 


24.3  CASE  (b).  SIMPLE  AND  BINARY  GAS 

The  procedure  is  essentially  the  same  as  before,  but  to  determine  the 
1st  approximation  we  now  get  the  equation 


(fix  w)  •  —  +  S  (f(o>,  f<o))  =  0 

dC 


(4.18) 


for  a  simple  gas,  and 


(C*  “  '  IT  *  pi,  s.e  (,io)’  f(*o)>  ■  0 


(4.19) 
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for  a  binary  gaa .  Thia  again  leada  to  Maxwellian  distribution  in  the  let  ap¬ 
proximation.  For  the  next  approximation  fO)  =  fW  f  0)  we  again  get 
the  expreaaiona  (4.7)  for  a  eimple  gaa  and  (4.  12)  for  a  binary  gaa.  Aa  op- 
poaed  to  the  previoua  caaea,  however,  A  ,  B  and  D  get  the  forme  (3.70) 
to  be  expected  for  a  rotationally  aymmetric  medium, 

A  =  <B (C)  •  C  ,  D  =  )r  (C)  •  C  ,  B  =  |p  (C)  :  C  °  C  ,  (4.20) 

where  « ,  J3  and  ft  are  all  rotationally  aymmetric  tenaora.  According  to 
Section  22  we  can  then  write  A  and  D  aa 


A  =  ^(CJC  +  «n(C)Cx  -°+  am(C)(Cx  «°)  x  *°  . 
D  =  /(CJC  +  rn(C)C  x  **°  +  rm(C)(Cxe°)  x 


(4.21) 


ov  o 

M 


while  B  ia  a  linear  combination  of  the  aix  independent,  divergency -free 
and  aymmetric  tenaora  that  may  be  found  by  combination  of  the  vectcra 

C  ,  C  x  «e°,  (  C  xa°)  x  e»°  .  These  tenaora  are  given  by  the  expreaaiona 
(3.72a).  The  coefficienta  a,  0  and  y  can  be  evaluated  in  the  same  way  as 
previously,  but  it  ia  evident  that  it  would  be  very  complicated  to  carry  out 
the  evaluations  in  practice.  The  procedure  ia,  of  course,  the  same  if  the 
number  of  components  of  the  gaa  ia  greater  than  two. 

Now  having  formally  determined  the  1st  and  2nd  approximations  to  the 
distribution  functions,  we  shall  uae  thia  in  the  following  section  to  find  the 
expressions  for  the  transport  vectors. 


25.  EVALUATION  OF  THE  TRANSPORT  VECTORS 


The  transport  vectors  or  fluxes  in  which  we  are  interested  are,  aa  be¬ 
fore,  diffusion,  current  density,  transport  of  momentum  and  the  heat  flow, 
i.e.  the  quantities 


n 


=  e  n 

«  a 


o 


n  m  C  °  C  , 

a  a  a  a 


n  ^m  C 
t  a  a  a 

(4.22) 


If  the  series  of  approximations  to  fa  converge  sufficiently  rapidly,  which 
ia  the  case  if  the  temperature  ia  sufficiently  small  and  the  density  is  suffi¬ 
ciently  large  (4.  1),  the  two  first  terms  in  the  series,  f^°)  +  f^  will  give 
a  good  approximation  to  fa  .  The  function  f^  =  f^°)  t  ^  is  given  by 
(4.  12).  The  contribution  to  (4.22)  from  the  1st  approximation  fj^  is  aero, 

as  all  these  quantities  are  aero  in  the  state  of  equilibrium.  The  contribu¬ 
tion  from  the  2nd  approximation  is 
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(a) 

JO 

9aC 

1 

II 

d  In  T  i  * 

dr  *k  *2 

(b) 

o(0 

pa  = 

-2«*a  : 

e 

(c) 

-  M 

a 

dlnT 

dr  «  12 

<d) 

II 

.  i'  . 

d  In  T  - 

dr  «  12 

where  the  coefficients  are  defined  by 

-  H  *  C«^(«O)d0fl  . 


r(°l 


-2*ia  =  JmeCttCa  ®a  Vdc.  . 
-»«  =  /imacJcaAa  4°)doa  * 

-o'  =  Je  C  A  f^do  , 
a  *  a  a  a  a  a 


(4.23) 


*.  =  /  m  C  D„  f^dc  , 

a  v  a  a  a  S  a 


•  i  Kc«wio)d<>.  • 

-•  *  J«JW(O)d0a  ' 


(4.24) 


For  A  ,  IB  and  D  we  further  substitute  the  quantities  (4.20).  It  is  then 
easily  shown  that  the  six  second-order  coefficient  tensors  included  in  the 
equations  (4.23  a-c-d)  have  precisely  the  expected  symmetry  properties, 
i.e.  they  are  rotationally  symmetrical. 

The  corresponding  circumstances  are  not  so  easily  shown  in  the  case 
of  equation  (4.23b),  but  evaluations  which  will  be  omitted  here,  show  that 
the  fourth -order  tensor  is  also  rotationally  symmetrical.  When  the  magne¬ 
tic  field  disappears,  A,  B  and  D  are  given  by  (4.  17)  and  all  coefficients 
(4.24)  will  be  isotropic  tensors. 

The  equations  (4.23)  which  are  derived  by  microscopic  theory  there¬ 
fore  are  completely  equivalent  with  the  corresponding  phenomenological  set 

(3.29-30).  The  choice  of  other  "forces"  (  4  ,  «  ,  d#  (4.15))  is  of 

course  unimportant,  as  they  are  linear  combinations  of  our  original  "forces" 

(AiUJL  ,  e  ,  ^Ji).  The  coefficients  m  ,  JA,  0  and  ut  in  Section  19  are 
0r  or 

therefore  given  linear  combinations  of  the  coefficients  (4.24)  (especially 
|i  s  (£  ).  They  can  therefore  be  evaluated  when  A,  IB  and  D  are  deter¬ 
mined  in  the  manner  sketched  in  the  previous  section. 

In  the  following  we  denote  the  respective  non -magnetic  and  magnetic 
cases  by  (a)  and  (b)  as  before. 


n 


1 
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25.1 

CASE  (a)  AND 

(b).  SIMPLE  GAS 

In  this  case 

d12 

=  0  and  (4.23) 

gets  the  form 

(a) 

-  (1) 

QC 

0 

(c) 

o 

II 

> 

(b) 

P(1>  - 

-  2*« i 

:  e  (d) 

(4.25) 

dr 


If  we  substitute  for  A  and  B  the  values  obtained  by  breaking  off  the  series 
of  Sonine  polynomials  (4.9)  after  the  first  term,  we  get,  according  to  Chap¬ 
man  and  Cowling  [4]  for  the  heat  flow, 


-(«)i 


: — 21 

1+W  T6 


/  d  In  T  , 

dr  'l 


WT 


l+i 


,2t2 


d  In  T 
d  r 


o  ,  ,  din  T 
x  td  +  (- 


dr 


y 


(4.26) 


Here  (x  ),  denotes  the  first  approximation  to  the  heat  conduction  coefficient, 
and  t  is  defined  as 


T 


Zm 
5  kp 


(3.27) 


For  solid,  spherical  molecules  this  is  identical  to  the  collision  time  for 
the  molecules,  apart  from  numerical  factors  [4]  .  Even  when  the  molecules 
are  electrically  charged,  as  in  our  case,  can  we  give  a  similar  interpreta¬ 
tion  to  t  .  The  factor  wt  is  therefore  essentially  the  ratio  between  the  col¬ 
lision  time  and  gyration  time  in  the  magnetic  field. 

Equation  (4.26)  means  that  the  heat  flow  coefficient  is  given  by  the 
matrix 


(>*>!  =  (*)] 


1+  wZ  xZ 


-  WT 

1  +  U)2  T2 


WT 

21 


1+  Ul  T 


2  2 
1+  U)  T 


(4.28) 


The  form  of  this  matrix  is  precisely  as  expected.  We  see,  for  instance, 
that  the  non -diagonal  and  anti -symmetric  elements  are  linear  in  the  magne¬ 
tic  field,  in  conformity  with  the  Onsager  relations  (3.76). 

We  see  that  the  factor  wt  is  a  measure  of  the  anisotropy  of  the  system. 
If  we  put  wt  =  0  ,  then  (4.28)  gives  tne  value  we  must  have  for  (ix)  j  in  the 
non -magnetic  case,  i.e.  (<x)j  =  (x)j  I. 

We  see  from  (4.26)  that  the  heat  flow  parallel  to  a  magnetic  field  will 
not  be  influenced  by  the  field.  When  wt  grows  from  zero  towards  unity, 

d  In  T 

the  flow  in  the  direction  of  -( — g-j —  )  will  gradually  decrease,  while  a 
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"Hall  flow"  is  a  set  up  in  the  direction  - 


din  T 

dr 


This  flow  is  at  a 


maximum  when  wt  =  ]  and  decreases  like  the  usual  transverse  current 
when  wt  grows.  For  very  strong  magnetic  fields  ( wt  »  1)  we  must  ex¬ 
pect,  according  to  (4.26),  that  the  heat  flow  across  the  field  becomes 
completely  subdued.  The  reason  for  this  is  of  course  that  the  effective  mean 
free  path  across  the  field  in  this  case  no  longer  is  X  =  c  •  t  but  the  Larmor 


radius 
lines . 


p  1 

=  —  =  —  .  The  usual  free  path  is  "curled  up"  around  the  field 

(0  U)T 


Chapman  and  Cowling  further  give,  without  proof,  the  first  approxima¬ 
tion  to  jp  with  the  remark  that  the  result  is  of  little  interest  because  the 
corresponding  phenomena  cannot  be  observed  in  metals  (!).  The  result  is 
given  on  a  component  form  which  is  not  easily  read,  but  closer  scrutiny 
shows  that  the  results  correspond  to  equation  (4.25)  with  the  coefficient 
tensors 
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The  1st  approximation  to  the  viscosity  coefficient  (p)j  is  given  by 


(4.29) 


<“>1  =  3  PT 


(4.30) 


where  *  is  defined  by  (4.27)  as  before. 

The  form  of  the  matrix  is  again  as  expected.  The  non -diagonal  and 
anti -symmetric  elements  are  linear  in  the  field.  The  Onsager  relations 
(3.76)  are  thus  once  more  satisfied.  If  we  write  wt  =  0  equation  (4.29) 
gives  the  result  obtained  for  2(|t),  in  the  non-magnetic  case  (Section  21). 
At  the  opposite  extreme,  wt-+  »  ,  all  elements  in  (4.29),  except  the  nine 
elements  in  the  upper  left-hand  corner,  approach  zero,  and  the  coefficient 
matrix  once  more  becomes  isotropic.  The  stress  tensor  now  becomes 


/i\  (i\  o  o  o 

P11  =  P22  =  '  (**  ( eli  +  e22^  =  *33  ' 

<i  ♦  j)  • 


P(I)  =  -2(n)  e 
p33  MP/j  *33 


p0>  =  0 
pij 


(4.31) 
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o 

If  i»  zero,  which  ia  the  case  if  the  gas  moves  along  the  magnetic  field 

with  the  velocity  wQ  =  wQ(  )  the  viscous  forces  will  completely  disappear. 

It  has  been  discussed  whether  or  not  the  derived  value  of  the  viscosity 
tensor  is  correct.  There  can  scarcely  be  any  doubt  that  the  result  is  correct, 
and  it  is  even  supported  by  other  works.  It  is  another  matter  if  one  can  sim¬ 
ply  assume  that  (4.29)  is  valid  when  wt  »  1.  This  appears  to  conflict  with 
the  assumptions  for  the  validity  of  the  method  of  solution  used  in  Section  24, 
and  corresponds  to  a  situation  where  the  term  (Cxs)  '  ^  in  Boltzmann's 
equation  is  dominant,  and  the  collision  term  must  be  regarded  as  a  pertur¬ 
bation  . 


25.2  CASE  (a)  AND  (b).  BINARY  GAS 
According  to  (4.23)  the  relative  diffusion  is  given  by 


Ci  -  c2 


r  .  ±}*jl  +  &  .  d  i 

alQ2  L  1  Or  1  12 J 


(4.32) 


as  the  relations  I  XV  =  Z  %  =  O  must  be  satisfied.  By  constant  tem- 

<*  “  a  * 

perature,  the  1st  approximation  to  (4.32)  is,  according  to  Chapman  and 
Cowling, 


(4.33) 


Here  (  A.j)j  is  the  1st  approximation  to  the  diffusion  coefficient,  and  w 
and  t  are  given  by  the  expressions 


1  , 

“  =  g  <*2  “I 


+  <Jl“2)’ 


*  * 


Thus  «  is  an  average  gyro  frequency  for  the  two  component  gases .  In 
a  similar  way  we  can  show  that  for  solid,  spherical  molecules  t  is  the 
average  collision  time  for  the  two  components . 

We  can  continue  to  give  the  same  interpretations  to  t  when  the  mole¬ 
cules  are  electrically  charged.  As  the  current  density  j*  is  proportional 
to  (4.33)  according  to  (3.8),  we  immediately  find  the  expression 


MT 


l+»2  T2 


d12x“ 


+  (<* 


(4.34) 


This  "Ohm's  law"  shows,  together  with  (4.33)  that  both  the  diffusion 
coefficient  tensor  and  the  conductivity  tensor  have  the  expected  form  (4.28). 

Chapman  and  Cowling  do  not  give  an  expression  for  the  heat  flow  and 
the  viscosity  tensor  in  the  magnetic  case,  but  we  must  expect  that  the  co¬ 
efficient  tensors  have  the  form  (4.28-29)  even  in  this  case. 
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The  expressions  found  in  the  non-magnetic  case  are  in  complete  accord 
with  (4.23)  for  an  isotropic  medium.  The  expressions  (4.33)  fcnd  (4.34)  re¬ 
duce  to  the  corresponding  expressions  for  the  non-magnetic  case  when  tst 
is  put  equal  to  zero. 


U 
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APPENDIX 


The  summational  invariants  $  of  the  Landau  collision  term  (2.  18)  are 
characterized  by  having  the  property 


o  , 

(A.  1) 


where 


,r->,  ■  d« 


n.  ( 


They  may  be  found  in  the  following  manner. 


2  2 

Using  equation  (2.  18)  we  get  by  partial  integration  (xap  =  2*Xeg  e^) 
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and  analogously 

d.fo  i  dfft  ,  d  Inf  dlnf' 

■  * !l *7  t**  •  u  ■  (S7  -15s-  ■  v.dc'de 

(A.  3) 


Summing  these  two  equations  we  get 
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(A.  4) 


It  is  at  once  evident  that  only  the  three  molecular  properties 


4>  =  1  (Const.),  me  ,  i 


me 


(A.  5) 


will  make  this  integral  and  n 
summational  invariants . 


*c4* 

dt 


zero.  Consequently,  these  are  the 


(We  also  notice  that  substitution  of  $  s  *  in  (A. 4)  at  once  gives 
us  the  proof  of  the  H -theorem.) 


m 


-  91  - 


ACKNOWLEDGEMENTS 


I  am  indebted  to  Dr.  E.  Jensen,  who  has  been  my  adviser  during  the 
work  on  this  paper,  and  to  Dr.  S.  Tjptta  for  some  clarifying  discussions. 

The  research  reported  in  this  document  has  been  sponsored  in  part 
by  the  Geophysics  Research  Directorate  of  the  Air  Force  Cambridge 
Research  Laboratories,  Air  Force  Research  Division,  UNITED  STATES 
AIR  FORCE,  under  Contract  No.  AF  61  (052)- 186  through  the  European 
Office  (EOARDC)  in  Brussels,  Belgium. 


M 


-  93  - 


REFERENCES 

[1]  Bogoliubov,  N.  N.  (1946)  -  Problems  of  Dynamical  Theory  in  Statisti¬ 

cal  Physics.  USSR. 

[2]  Brittin,  W.  E.  (1957)  -Phys.  Rev.  106:843. 

[3]  Holter,  6.  (1959)  -  On  the  Propagation  of  Transverse  Electromagnetic 

Waves  in  a  Plasma  in  a  Magnetic  Field.  Report  No.  7.  Inst,  of 
Theor.  Astrophys  .  ,  Oslo. 

[4]  Chapman,  S.  and  T.  G.  Cowling  (1939)  -  The  Mathematical  Theory  of 

Non-Uniform  Cases,  Cambridge. 

[5]  Landau,  L.  D.  (1937)  -  JETP.  7:203.  USSR. 

[6]  Bohm,  D.  and  D.  Pines  (1952)  -  Phys.  Rev.  85:338. 

[7]  Kadomtsev,  B.  B.  (1958)  -  JETP.  6:1.  USSR. 

[8]  Green,  H.  S.  (1959)  -  Phys.  of  Fluids.  2:4. 

[9]  Enoch,  J.  (I960)  -  Phys.  of  Fluids.  3:3. 

[10]  Fowler,  R.  H.  (1936)  -  Statistical  Mechanics ,  Cambridge. 

[11]  Boltzmann,  L.  (1872)  -  Wien.  Sitz.  66:275. 

[12]  Lorentz,  H.  A.  (1887)  -  Wien.  Sitz.  95:127. 

[13]  Jeans,  J.  (1916)  -  The  Dynamical  Theory  of  Gases,  Cambridge. 

[14]  Braginskij,  S.  I.  (1958)  -  JETP.  6:33.  USSR. 

[15]  Spitzer,  L.  jr.  (1956)  -  Physics  of  Fully  Ionized  Gases ,  New  York. 

[16]  Chandrasekhar,  S.  (1943)  -  Astrophys.  J.  97:255. 

[17]  Cohen,  R.  S.  ,  L.  Spitzer  jr.,  P.  McR.  Routly  (1950)  -  Phys.  Rev. 

80:230. 

[18]  Schluter,  A.  (1957)  -  Zs.  f.  Naturforsch.  12a:822. 

[19]  Schatzmann,  E.  (1959)  -  Nuovo  Cimento,  Suppl.  13. 

[20]  Spain,  B.  (1953)  -  Tensor  Calculus,  Edinburgh. 

[21]  de  Groot,  S.  R.  (1951)  -  Thermodynamics  of  Irreversible  Processes, 

Amsterdam. 

[22]  Milne -Thomson,  L.  M.  (1955)  -  Theoretical  Hydrodynamics,  London. 

[23]  Chew,  G.  F.  ,  M.  L.  Goldberger  and  F.  E.  Low  (1956)  -  Proc.  Roy. 

Soc.  A.  236:112. 


n 


